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ABSTRACT 


The Ising model of a ferromagnet has been 
investigated in detail on a triple of related lattices 
which heretofore have not been studied. These lattices, 
which we call the hydrogen peroxide, hyperkagome, and 
hypertriangular have coordination numbers of three, four, 
and six respectively, and may be looked upon as three 
dimensional analogues of the more familiar two dimensional 


lattices, honeycomb, kagome, and triangular. 


High and low temperature exact series expan- 
Sions have been derived for the partition function, 
energy, specific heat, and initial susceptibility of the 
hydrogen peroxide and hypertriangular lattices. The exact 
series expansion of the spontaneous magnetization has also 


been derived for the hydrogen peroxide lattice. 


The star - triangle and decoration transforma- 
tions, which relate the Ising partition functions of the 
hydrogen peroxide, hyperkagome, and hypertriangular lattices 
te each other, have been utilized to obtain the properties 
of the hyperkagome lattice from those of the hydrogen 
peroxide lattice. For a given pair of the triple of related 
lattices we have utilized these transformations to obtain 


expressions relating the specif heat amplitudes and the 
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susceptibility amplitudes. 


A problem which is very similar to the calculation 
of the initial susceptibility of an Ising ferromagnet, 
namely the self-avoiding walk problem, has also been 
studied for the hydrogen peroxide and hypertriangular 
lattices. Exact series expansions of the self-avoiding 


walk generating function have been derived for both lattices. 


A numerical analysis utilizing the ratio and 
Pade approximant methods has been performed on the series 
expansions to obtain estimates of the critical constants 
and critical indices of the Ising model. Estimates for the 
Critical index and the attrition parameter of the self- 


avoiding walk problem have been obtained in a similar manner. 


An empirical formula giving the critical point 
of an Ising ferromagnet as a function of coordination 


number has also been found. 
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CHAPTER I 


INTRODUCTION 


A. Critical Point Phenomena 


The study of critical phenomena was initiated 
a century sue with the measurements by Andrews (1869) near 
the critical point of carbon dioxide. Subsequent 
experiments have shown that a number of physical systems of 
a diverse nature, fluids, ferromagnets, antiferromagnets, 
binary liquids, and binary alloys, all have well defined 
critical points and all show very similar behaviour in the 
vicinity of the critical point. One could also mention 
critical point phenomena associated with superfluids, 
superconductors, ferroelectrics and antiferroelectrics. 
The origin of critical phenomena in these systems is how- 
ever of’a.different nature from those systems we have 


mentioned above so we will say no more about them. 


The major part of this thesis is concerned with 
the theory of. critical phenomena described in terms of an 
Ising model of a ferromagnet. Hence it is necessary that 
we first give a brief review of the experimental facts as 
we know them. An excellent review of the methods and 


results of experimental investigations of critical point 
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phenomena has been given by Heller (1967). We will limit - 
most of our discussion of experimental results to critical 
phenomena associated with gas-liquid and magnetic systems ~ 
as they have been subject to more extensive investigations 


than other systems. 


Our discussion of the experimental and theoretical 
aspects of critical phenomena will concentrate on the 
manner in which various thermodynamic functions behave as 


the temperature approaches its critical value. 


Experience has shown that in the neighbourhood 
of the critical point certain thermodynamic functions behave 


as 


x * Pp a 
W 2 A, (x X,) eS Ke 5 (1-1) 


where W is the thermodynamic function of interest and 
x is the independent thermodynamic variable. The exponent 
p is called the critical exponent or critical index of 
Wheand« Ay is called the amplitude. In case W diverges 
at the critical point we will write 

W = Bea dren. bh, xast ix ‘ (1-2) 
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will also occur. In each, case p is called the critical 
index of the thermodynamic function W . It has become. 
established notation to specify by Greek letters the 
critical indices of certaii thermodynamic functions, e.g. 
@ for specific heat,- 6: for-magnetization,. y ‘for 
susceptibility etc. There is one other form of critical 


point behaviour which often occurs in the critical region, 


WwW = A 1n(x-x,) +B, x> oo (1-4) 


A logarithmic divergence such as (1-4) can be described by 


a eritieal index p = 0 


The phase diagram for a typical single component 
system is depicted in Fig. 1.1. The point labelled xX at 
the end of/the vapour pressure curve, beyond which it is 
impossible to make any distinction between the liquid and 
gas phases is called the critical point (PSP oy T=T.)- A 
transition from the gas to the liquid phase along curve 


AB:-(Fig. 1.1) is continuous and does not give rise to any 
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Fig. 1.1. The phase diagram for a simple 
substance exhibiting gas, liquid, and solid 
phases. The point labelled xX is the 


eritical point (T: >= Toe P = P) 
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anomalous behaviour in the derivatives of the free energy,” 
while the transition from the gas to the liquid phase along 


curve CD (Fig. 1.1) is discontinuous. 


A more convenient description of the critical 
region of a fluid and the phenomena associated with it can 
be given in terms of the shape of the isotherms in the 
region ef the. critical. point (Fig. 1.2). The pressure is 
a monotonically increasing function of the density for 
temperatures greater than the critical temperature To 4 
and its derivatives are finite and continuous. As T is 


Cc 


approached from above (T > To ) ae (OP/90),, decreases and 
Sa 
has’a minimum for some value of pop until at Te=T and 


(©; 
pee Par (OP/OP) , = 0. As the temperature is lowered 
below To the fluid separates into two phases, gas and 
dciquid;, and (9P/39) 5 vanishes over the interval 


Pa < OS Oy where is the maximum density of. the gas 


We: 


phase and is the minimum density of the liquid phase. 


we 
The shape of the isotherms in the critical region 


can be characterized as follows. For simple gases 


(Ar; Xe, CO , ect.) the top of. the coexistence curve can be 
2 


described by the form 


A (T-T) ’ T+ T ° (I=5) 
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Fig. 1.2. The isotherms of a simple fluid. 
The coexistence curve is indicated by the 
dashed line and the critical point is at 
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An analysis of the data from the classic experiment on 


Xe (Weinberger and Schneider 1952, Fisher 1964(a)) yields 
8-= 0.345 +0.015 - 


The shape of the critical isotherm near the 
critical point is described by 
i 


Bee gan ere « 


, ee 8) pal (1-6) 
Widom and Rice (1955) analyzed experimental data for a 


number of simple gases. They found 


The appropriate quantity to describe the shape of 
‘the critical isotherms above To is the isothermal 


compressibility, defined by 


Lev 1 ,dp 
K, = - Slyp) = S(x) (1-7) 
TS Viv or nm oe dp T 


which according to our above statements approaches infinity 


aS. Te T. . The behaviour of the isothermal compress- 
+ 
ibility in-the critical region.can be described by 
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+ 
The experimental evidence for CO. (Heller 1967) suggests 
Vel 35 £0.25 


The behaviour of the compressibility along the coexistence 
curve as the critical point is approached from below can 
be similarly described, 


a 


Re (Palo Sete.) oe Chai) oo (1-9) 


An analysis of the data for co. (Heller 1967) indicates 


Recent experiments (Voronel' et al. 1965, 
Moldover 1966) have revealed a very striking anomaly 
associated with the specific heat at constant volume as 
one approaches the. critical point from above and below 
along the critical isochor (D=—,)- The experimental 


evidence (Kadanoff et al. 1967) suggests 
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where 


A value of a (or a4~) equal to zero would indicate a 


logarithmic divergence. 


Guggenheim (1945) has shown that most simple 
gases*obey a law of corresponding states. This suggests 
that the values of the critical indices listed above 


Should be valid for all such gases. 


The phase diagram of an idealized ferromagnet 
is illustrated in Fig. 1.3. There are two phases, an 
ordered phase represented by the heavy line on the T axis, 
and -a disordered phase represented by the remainder of the 
diagram. The ordered phase, for which there exists a 
spontaneous magnetization, occurs only for H = 0 and 


T aap . The point labelled X(T=T) is the ferromagnetic 
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Fig. 1.3. The phase diagram of a 
ferromagnet. The point X(T=T, H=0) is 


the ferromagnetic Curie point. 
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The critical point phenomena of a ferromagnet 
can also be explained in terms of the shape of the isotherms 
(Fig. 1.4). For a given direction of the magnetic field 
H and, temperatures greater than the critical temperature 
TO the magnetization is a monotonically increasing function 
of the field and has finite and continous derivatives for > 
gali,-valuesof:.the fielid. For Tos To the isotherms 


show a discontinuity at H = 0 , the magnetization tending 


to the limits I. (0) and I-(0) as’ H tends-to zero. 


Near the critical point the spontaneous magnetiza- 


tion behaves as 


Cc 


t . (p-0)°, 1c Se Oe oe eae (fea) 


For the insulating ferromagnet Eus (Helier and 


Benedeck 1965) one finds 
Ss O833 £0.02 


The shape of the critical isotherm can be described 


by 
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Fig. 1.4. The isotherms of an idealized 


ferromagnet. 
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ena eal ele Ste 19.0 (1-12) 


Kouvel and Fisher (1964) analyzed some experimental data 


for nickel: and found 


oe GIES) 


increases gradually as the temperature is lowered until at 
T = To and H= 0 it becomes infinite (Fig. 1.5). The 
behaviour of the initial susceptibility in the region of 


‘the critical point is described by 


ie me, , = 0. Siar 
pil a aia A TH 0 (1-14) 


A very accurate experiment on Fe (Noakes et al. 1966) 


indicates 


y = 1.333 + 0.015 
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Fig. 1.5. The spontaneous magnetization 


and initial susceptibility of a ferromagnet. 
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One would expect similar behaviour for the initial susceptibility 
below T 

Cc 
eek 


a Cc a 


eee oer ye H=0, (1-15) 


At the present time there are no experimental results which 


@ 


enable one to predict a value for y 


The critical point behaviour of the specific heat 


Cu at constant magnetic field is described by 
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Cy (psa ie : Dat ge 
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Experiments performed by Teaney (1966) show the ferromagnets 
EuS and Eu0 both have magnetic specific heat anomalies 

of the lambda type. The experimental results are not 
sufficiently accurate to allow one to fix a definite value 
for o- and oa* . All that can-be said is that the results 
are- consistent with a logarithmic behaviour or a small non 


zero value- for the. critical index. 


In simple antiferromagnets the quantity analogous 


‘to the spontaneous magnetization of a ferromagnet is the 
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sublattice magnetization. Experiments performed on the 


antiferromagnet MnF, by Heller and Benedeck (1962) give 
8B = 0.335 + 0.005. 


Teaney (1965) has made accurate measurements of 
the specific heat of the antiferromagnetic MnF, . The 
behaviour of the specific heat in the critical region 
(Kadanoff et al. 1967) is consistent with 

a 


Go (tate) n+ nate 


HO, (1-17) 
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In order to have a critical isotherm and 
susceptibility of an antiferromagnet analogous to those 


of a ferromagnet it is necessary to place the antiferro- 
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magnet in-a staggered field, a field which points in 
opposite directions for neighbouring spins. There is no 
known way of realizing such a field experimentally, hence 
there is no direct experimental analogue between the 
critical isotherm and susceptibility of a ferromagnet and 


antiferromagnet. 


Many binary fluid mixtures exhibit a critical 
temperature below which the fluid separates into two 
phases. One can define a coexistence curve analogous to 
that of a gas-liquid system. Some recent accurate measure- 
ments on mixtures of C Cl and n- Ca (Thompson and 
Rice 1964) indicate that the index 8 describing the 
shape of the coexistence curve in the critical region is 


given by 


One other type of transition which we shall 
mention here is the order-disorder transition in a binary 
alloy. A binary alloy in which like atoms attract unlike 
atoms is analogous to an antiferromagnet. The quantity 
analogous to the sublattice magnetization is the order 
parameter. A very accurate experiment has been performed 


on very pure beta brass (Als-Neilsen and Deitrich 1967). 
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An analysis of the experimental data suggest 


B = 0.305 + 0.005 


The behaviour of the diffuse scattering at the superlattice 
peak was also investigated. This provides information. on 
the analogue of. the staggered susceptibility. The results 


showed 


In light-of- the results:ffor § and y found 
from our work on the Ising model the above results for 
beta brass are perhaps the most significant of all the 


experimental results we have quoted. 


The experimental evidence we have cited is 
summarized in Table 1.1. The values of the critical 
indices suggest that phenomena associated with critical 
points are very general and display only a very weak 
dependence: on a particular type of system, i.e. magnetic, 
gas-liquid, etc. Thus it should seem quite reasonable 
to seek- expressions relating the critical indices to one 
another. Using rigorous thermodynamic arguments, the 


relation between the specific heat at constant field Cy 
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and the specific: heat at constant magnetjzation Cur and: the 
fact that Cu is non-negative, Rushbrooke (1963) derived 
the inequality 


2 


Gi two gece. yoorg 2 te (1-18) 


Griffiths (1965) gave a rigorous derivation, based on the 


convexity properties of the free energy, of the inequality 
OY  F0B41+G]oatoor (12h) 


A-non-rigorous. theory has. recently been 
developed -(Widom 1965, Kadanoff 1966, Kadanoff et al. 1967), 
which enables all of the critical indices to be. expressed 
in terms of two fundamental indices. The theory (scaling 
law theory) predicts that the inequalities in (1-18) and 
(1-19) are replaced by equalities and that corresponding 
indices’ above and below the critical point are equal, i.e. 
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The experimental values for the set of critical 
indices describing the critical point behaviour of real 
systems are not known to sufficient accuracy to enable. one 
to draw a definite conclusion on the validity of. the 
scaling law hypothesis which leads to (1-20). The critical 
indices of the mean field theory and also those of the 
two dimensional Ising model satisfy (1-20). However, 
numerical investigations indicate’(1-20) does not hold for 
the three--dimensional Ising model (see Table 1.2). We 
shall have: more to say about equations (1-18), (1-19), and 


(1-20):-in-. Chapter VIII. 


Several attempts. have been made towards a 
theoretical description of the experimental facts as out= 
lined above. -The first of these were of a phenomenological 
nature, the Van der Waals theory of the gas-liquid critical 
point and the Weiss theory of ferromagnetism, which are now 
commonly known as mean field theories. The predictions of 
the-mean field theories disagree both qualitatively and 


quantitatively with the results of experiment. 


An explanation based on a statistical mechanical 
approach did not come until much later with the advent of 


the Ising (1925) and Heisenberg (1928) models of ferro- 


magnetism. 
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Theoretical investigations of the Ising and 
Heisenberg models have shown these two models, or modifica= 
tions thereof, are capable of explaining much of the 
observed critical point behaviour in real systems. Our 
work concentrates on the Ising model of a ferromagnet, so 
we will review in the next section the development of 


the Ising model from 1925 to the present. 
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B. The Ising Model 


What has now come to be known as the Ising: model 
was first proposed as a model of ferromagnetism in 1920 to 
Ernst Ising by his supervisor W. Lenz (Brush 1967). The 
Ising model represented the first attempt to give a: 
statistical mechanical interpretation to the phenomena of 


ferromagnetism. 


The classical Ising model of a ferromagnet 
associates with each site of a crystal lattice a spin 
(magnetic moment m) which can point in two directions, "up" 
and "down", relative to an external magnetic field H. The 
two states of the ith spin are usually described by the 
variable is +1. for "spin-up" and a, = =i “for, "Spin 
down". In addition to the energy of- interaction (-mHo . ) 
of the magnetic moment with the magnetic field, there is 
also assumed to exist an interaction between neighbouring 
spins, -dJ if the two spins point in the same direction and 
J if they point in opposite directions, where J is a 
positive constant. The energy of the system is given by 


E.e--J }. o;9, -mH joo, (1-21) 


ae al 4. 


where <i,j> indicates the first sum is over nearest 


neighbours only.» Equation (1-21) can be taken as the 
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definition of the classical Ising model. Having defined the 
energy the next step is-to calculate (in principal: at: least) 
the partition function from which the thermodynamic 
quantities of interest can be obtained in a straight forward 
manner. The greater part of this thesis is concerned with 
just such a calculation for a triple of related three 
dimensional lattices on which the Ising model (1.21) has not 


previously been studied. 


Interest in the Ising model is primarily of a 
pedagogical nature as it is one of the few many body 
problems that has proven exactly soluble, others being a 
class of two dimensional problems (the ice problem, KDP 
model of a ferroelectric, and F model of an antiferro- 
electric) ,which have been solved by Lieb (1967(a), (b),.(e)), 
and the one’ dimensional Van der Waals gas,which has been 
solved exactly (Kac et al. 1963). Despite the fact the 
expression for the energy (1-21) is extremely simple, 
studies of the Ising model have shown that the simple nearest 
neighbour interaction between spins is the essential 
ingredient in explaining the existence and behaviour of 


ferromagnetic systems near their critical points. 


Use of the Ising model is not restricted to a 
study of ferromagnets only. The model can serve as an 


equally: valid description of an antiferromagnet upon 
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replacing J by -J3J , as. a model for liquid-gas° condensa- 
tion; the lattice gas (Yang and Lee 1952); and as a model 


for binary alloys and liquids as well. 


Though as we have pointed out in the previous: 
paragraph the great pedagogical interest in the Ising 
modei stems from the fact it serves as a prototype for many. 
different types of cooperative phenomena, there is experi- 
mental justification for studying such a model as well. For 
those insulating ferromagnets and antiferromagnets exhibiting 
a large anisotropy in the g factor (dy > gy) the Ising 
model provides a more valid description than the Heisenberg 
model which applies only to isotropic systems. Recent 
experimental evidence (Als-Neilsen and Deitrich 1967) 
suggests the Ising model can provide a valid description of 


a binary alloy. 


Ising (1925) was able to solve the model exactly 
in one, dimension for all temperatures and all values of the 
magnetic field by a straight. forward combinatorial approach. 
The one dimensional model did not exhibit a phase transition 
for T° > 0 . However, at T= 0, the susceptibility becomes 
exponentially infinite, indicating the onset of a spontaneous 
magnetization. In this sense the one dimensional model can 
be-said to have-a critical point at T= 0-. As a result 


of his. studies, Ising erroneously concluded that. the two 
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dimensional model would also fail to exhibit ferromagnetic 
properties. Ising had failed to appreciate the importante 


of correlations between spins in two or more directions. 


The first convincing argument. that the two: 
dimensional model would exhibit ferromagnetic properties 
was given by Peierls (1936). The first exact results for 
the two dimensional Ising model were obtained by Krammers 
and Wannier (1941) who discovered a certain symmetry ~ 
property between the high and low temperature partition 
functions for the square lattice which enabled them to locate 
the position of the critical point exactly if.it existed. 
An exact closed form expression for the zero field parti- 
tion function for the rectangular lattice was first given 
in the classic paper of Onsager (1944), where he showed that 
the specific heat was logarithmically infinite on both sides 


et the ‘critical point. 


Onsager (see Wannier, 1945) was also able to shed 
light’ on. the reason for the symmetry property first 
observed by Krammers and Wannier. Onsager pointed out that 
the observed symmetry was unique to two dimensional lattices 
and could be given a topological.explanation based on the 
fact that two dimensional lattices are planar and therefore 
dual lattices (Ore 1963) can be constructed from them. The 


square lattice is self dual, the reason that Krammers and 
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Wannier were able to locate the critical point exactly. 


In 1948 Onsager (see Onsager 1949) enunciated-his 
famous result for the spontaneous magnetization of the 
rectangular lattice 


478 
eS HP (sinh 2K sinh a , (1-22) 


where 


(1-23) 


A 
Ii 
Gq 
ke 
= 
na 
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The derivation of this result has never been published by 
Onsager,vonly quoted in the above reference. The first 
published derivation of (1-22) was given by Yang (1952). 
Following the work of Onsager and Yang exact expressions 
for the zero field partition function and spontaneous 
magnetization for several other two dimensional lattices 
were obtained. The results, however, were of the same 
general form as those forthe square lattice, the critical 
indices in fact were found to be the same for all two 


dimensional lattices. 
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To this date no exact expression has been obtained 
for the partition function of the two dimensional Ising 
model in a non zero field, and hence no exact expression 
has been obtained for the initial susceptibility. “By 
examining the relationship between the correlations for 
pairs of spins and the initial susceptibility Fisher (1959(a)). 
was able to show in a semirigorous manner that near’ the 


critical point the initial susceptibility behaved as 


le SOL i a GR PN ear CO (1-24) 


The three dimensional Ising model has so far 
resisted: all efforts to obtain an exact solution, in.fact 
an exact result of any kind has. yet to be found. The only. 
rigorous results obtained so far for the three dimensional 
Ising model are those of Griffiths (1964-67) which show 
that. the three dimensional Ising model does indeed undergo 


a transition. 


In the absence of any exact analytical results 
for the three dimensional Ising model, most efforts at 
elucidating the behaviour of the thermodynamic: functions in 
the critical region have been limited to two alternative 
_methods, improved closed form approximations and exact 


series expansions. 
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In order to construct. a closed form approximation 
it is necessary to. neglect some of the statistical mechanical 
details of the problem so that the partition function can. 
be summed. The disadvantage of this approach is that the 
details» which one is forced to neglect are usually all 
important. Hence one simply reproduces the already well 
‘known mean field results which predict the wrong type of 
behaviour, both qualitatively and quantitatively, in the 


critical region. 


Our approach to the Ising model in this thesis 
is that of exact series expansions, the coefficients of the 
series being known exactly as far as they go. The method 
involves an exact: calculation of the leading terms (usually 
the first ten to twenty terms) of the power series expansion 
of the thermodynamic quantity of. interest. One then. attempts 
to estimate the critical parameters by employing suitable 
numerical techniques. The method has been tested on the 
two dimensional model and it has proven capable of giving 
excellent*estimates of the critical parameters, many of 
which are known:exactly. The impetus for this method was 
provided by Domb, Sykes, Fisher and their colleagues at 
Kings College,. London. As a result.of their efforts and 
those. of. others who have utilized this method a considerable 


amount of numerical information has been compiled on the 
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‘behaviour of the three dimensional Ising model in the région 


of the critical point (see Table 1.2). 


For example, in the region just above’ the: critical 


point it is now fairly well established that the initial 


ferromagnetic susceptibility behaves as 


= = fs = 
Ket t/t ke 1 HeO 


(1-25) 


There is also strong evidence to suggest that the 


specific heat: singularity above ZT, is not. logarithmic. 


The work of Sykes, Martin, and Hunter (1967) suggest that 


near To the specific heat behaves as 


ane. 
Cy = A, (1-T./T) ’ 


(1-26) 


On the low temperature side of To the values 


of the critical indices are much less certain. Probably 


the best previous estimates of the low temperature 


critical indices have been given by Baker and Gaunt (1967). 
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The primary purpose of this research has been to obtain: 
improved estimates of the low temperature critical indices. 
Such information is not only desirable from the standpoint 
of making a comparison with experiment but it can also. serve 
as a guide to the validity of the scaling law hypothesis 


(Widom 1965, Kadanoff 1966, Kadanoff et al. 1967). 
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€. Scope of. this. Thesis 


The structure of the three dimensional lattices 
on which we have applied the Ising model is discussed in 
Chapter II. The derivation of high temperature series 
expansions for the zero field partition function and initial 
susceptibility is outlined in Chapter III. The derivation 
of the low temperature series expansion for the partition 
function is discussed in Chapter IV. In Chapter V we 
review some well known Ising model transformations which 
are applicable to the lattices we have worked with. The 
self-avoiding walk problem, its relationship to the Ising 
model and the derivation of the corresponding series 
expansionsys is discussed in Chapter VI. Numerical methods 
utilized to obtain estimates of critical parameters are 
outlined in Chapter VII. A review and discussion of the 


results of this research project appear in Chapter VIII. 
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CHAPTER, LE 


THE HYDROGEN PEROXIDE, HYPERTRIANGULAR, AND 


HYPERKAGOME LATTICES 


Prior to the work reported on in this thesis, 
most investigations of the three dimensional Ising model 
have been confined to five regular lattices, the face 
centered cubic, body centered cubic, simple cubic, crystobalite, 
and diamond,having coordination numbers 12, 8, pe: and 
4 respectively. By a regular lattice we mean one which can 
be superimposed on itself by a combination of translation 
and improper rotation in such a way that any point can be 
superimposed on any other point and any bond can be superimposed 
On any other bond. In other words, in a regular lattice all 


points and all bonds are equivalent. 


The lattice of primary importance in this investi- 
gation is an idealized hydrogen peroxide lattice (Fig. 2.1). 
This lattice was first reported on in the literature by 
Heesch and Laves (1933) and a detailed description of this 
and other similar lattices appears in Wells (1954). This 
lattice belongs to the cubic system and has the space group 
0°(14 32) (Henry and Lonsdale 1952). The lattice is not 
Lae eee and so can occur in "left-handed" and "right- 


handed"*forms. The lattice consists of four interpenetrating 
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body centered cubic lattices (Fig. 2.2) and so has four atoms 
per unit cell. A coordinate system can be chosen such that 
the atoms in a cubic unit cell have coordinates (0,0,0); 

ey 2/450); (1/21/4174): (374,0;1/4); (0,374,374); 

4 243/5)5) (172,172,172) 3 (3/4,3/4,1/2) (Fig. 2.3): 
Idealized crystalline hydrogen peroxide is the only known 
substance to have such a structure (Wells, private communi- 
cation), hence the name of the lattice. The hydrogen per- 
oxide lattice may be constructed by picking out one quarter 


of the lattice points of a face centered cubic lattice. 


Real crystalline hydrogen peroxide differs from 
our ideal model in that it belongs to the tetragonal system 
(Abrahams et al. 1951) instead of the cubic system. The 
lattice parameters are a = 4,06 A and b = 8.00 A. The 
corresponding space group is Di (P4,2,2) (Henry and Lonsdale 
1952), enatiomorphous with DB(P,32,2). The oxygen-oxygen 
bond distance is 1.49 A, the distance between oxygen atoms 
across the hydrogen bond being 2.78 A. The angle made by 
the hydrogen bond with the oxygen-oxygen bond is 97° and 
the azimuthal angle between planes containing the hydrogen 
bonds and the oxygen-oxygen bond is 94° (Fig. 2.4). The 
hydrogen bonds form a four fold infinite helix in a plane 


perpendicular to the c-axis (Fig. 2.5). On the basis of 
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hydrogen peroxide lattice. 
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the structure outlined above, it was concluded (Abrahams 
aera i952) that crystalline hydrogen peroxide has no 
measurable entropy at absolute zero. For those interested. 
in a more detailed description of hydrogen: peroxide, from 
its discovery to the present, one may refer to the lengthy 


treatise by Schumb, Satterfield, and Wentworth (1955). 


In. studies of the three dimensional Ising model, 
the hydrogen peroxide lattice occupies a unique position. 
It is a lattice with the lowest possible coordination number 
(q=3) and it can be looked upon as the three dimensional 
analogue of the two dimensional honeycomb lattice. Further- 
more, from the hydrogen peroxide lattice it is possible to 
construct two other three dimensional lattices, the hyper- 
triangular and hyperkagomé lattices, having coordination 


numbers six and four respectively. 


The hypertriangular lattice can be formed from 
the hydrogen peroxide lattice by first bi-colouring the 
Sites of the hydrogen peroxide lattice (Fig. 2.6) followed 
by a removal of all sites belonging to one of the two given 
colours: The remaining sites form what we call the hyper- 
triangular lattice (Fig. 2.7). The hypertriangular lattice 
serves as a three dimensional analogue of the two dimensional 


triangular lattice. 
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The hyperkagomé lattice can be constructed from 
the hydregen peroxide lattice’in the following manner. At 
the mid peint of each bond of the hydrogen: peroxide lattice 
we place-a>new lattice site,° thus forming a decorated 
hydrogen peroxide lattice’ (Fig. 2.8). These’ new lattice 
Sites are then connected together and the original hydrogen 
peroxide sites are removed. The remaining lattice of 
coordination number four is called the hyperkagomé lattice 
(Fig. 2.9)... This lattice can be looked upon as being a 
three dimensional analogue of the two dimensional kagomé: 


lattice. 


The hydrogen peroxide lattice is a loose packed 
lattice whereas both the hypertriangular and hyperkagome 
lattices are close packed. By a loose packed lattice we 
mean a lattice in which there are no circuits (loops) made 
up of an odd number of bonds (edges). A close packed 
lattice has loops with both an even and an odd number of 


edges. 


One may be inclined to ask the following question, 
namely, why study the hydrogen peroxide lattice since it 
would- seem to be of little direct experimental interest? 

The hydrogen peroxide lattice, as we have already pointed 
out; occupies. a. unique, position in studies of the three 


dimensional Ising model. Because of its low coordination 
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number, the coefficient of the low temperature series 
expansions for the Ising model on the hy drogen-peroxide 

lattice are all of the same sign, a fact which tends to make 

| the lattice more suitable for study than lattices of higher 
coordination number. A notable exception of course being 

the diamond lattice with a coordination number of four (Essam 
and Sykes 1963), The low Pes een number of the hydrogen 
peroxide levthce also gives rise to a considerable simplifica- 
tion in the configurational counting problem which arises 

in the derivation of high temperature series expansions for 

the Ising model. Another reason for studying such a lattice 

is that the hydrogen peroxide, hypertriangular, and hyperkagome 
lattices are all related by certain well known transformations, 
which are described in Chapter V. These transformations 

serve as a valuable check as to the correctness of the 
underlying counting problem and they also provide an independent 
means of checking the results of the numerical analysis of 


the series expansions. 


The behaviour of the spontaneous magnetization, initial 
susceptibility, and zero field specific heat can all be 
interpreted in terms of their corresponding correlation func- 
tions (see Kadanoff et al. 1967). For example, the initial 
susceptibility can be formulated in terms of a sum over spin- 


spin correlation functions. As the critical point is approached, 
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the range of correlation grows, becoming infinite at the 
critical point, with the result that the susceptibility becomes 
infinite at the critical point. This diveregence in the range 
of correlation is the precise cause of this and all other 
singularities in the thermodynamic derivatives. Since one 
expects the long range correlations to be independent of all 
but the grossest features of the interaction, a change in the 
behaviour of the long range correlations from one lattice to. 
another of the same dimensionality is not expected. In fact 
the solutions for the two dimensional Ising model bear this 
out, the correlations for the square and triangular lattices 
being basically the same. In three dimensions as far as one 
can tell, the critical indices are also lattice independent. 
Hence, our numerical estimates of the critical indices for 
the Ising model on the hydrogen peroxide lattice are also 


expected to be valid for other three dimensional lattices. 
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CHAPTER III 


HIGH TEMPERATURE SERIES EXPANSIONS 


A. The Zero Field Partition Function 


In this chapter we briefly review some of the 
standard techniques used to derive high temperature series 
expansions. A thorough review of the various methods of 
deriving high temperature series expansions for the Ising 
model has been given by Domb (1960). It is not the 
intention of the author to give a lengthy review of the 
concepts of graph theory, which are widely used in deriving 
the expansions we are seeking. The reader unfamiliar with 
such concepts should refer to Domb (1960) or Sykes, Essam, 


Heap, and Hiley (1966). 


The energy of the spin one-half Ising model with 
nearest neighbour interactions only can be written in the 


form (c.f. Chap. I, eq'n. (1-21)), 


E=-JI oj o, -™ H z Os (3-1) 
<i,j> i 

where J is the interaction energy between neighbouring sites, 

m is the magnetic moment, and H is the external magnetic 


field. The o variables take the values + 1 according to 
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whether the magnetic moment is parallel or antiparallel to 


the magnetic field. 


The thermodynamics of the Ising model are computed 
from the partition function, which for a regular lattice of 
N sites takes the form 

Zy (8,H) = & exp B[J2ojo4+ mHZOi), (3-2) 

OyHt1,°o°° yet] 
where 8 = 1/kT and the outermost sum is over the 2N possible 
values of the o; for the N lattice sites. Since the o 


variables commute (3-2) can be written as a product 


Z =z > I Ko;.o4) Il ex mHO 3-3) 
nN (8/H) = E51, exp (Koyo) [ exp (6mio,), 


O,=+1, ie On =td 


where K = 8J/kT. The T4945 satisfy the following relations, 


ae soa)” = toe), © aie mete 
(ojo5)" = (oy95)° = (9595) a 
(3-4) 
I A 5 3 = ° ° -) S eee 
(0,05) = (0405) (5405) 
Hence we can write 
exp (Kojo4;) = cosh K + 040; sinh K. (3-5) 


Equation (3-5) is the well known Van der Waerden identity 
(Van der Waerden 1941). The first product in (3-3) can now 


be written as 
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(cosh K + ojo4 sinh K) = I cosh K (1 + 0305 tanh K) 
<i7,3? <i,j> 


gN/2 
= (coshK) , E[1+(ojo3) tanh K 
Fe ecNa 


+ (6495) (o,0,) tanh? K+---], 
(3-6) 


where qN/2 is the number of nearest neighbour bonds in the 
lattice. Setting H=0, the second factor in (3-3) becomes 


equal to unity and the zero field partition function is 


given by 
qN/2 
Zy (8,0) = (cosh K) ‘oie? {1 + (0405) tanh K 
at a 
Oy=1,°°° ,oy=tl (3-7) 


+ (9595) (0,04) tanh? K + ---], 


We can now give a graph theoretical interpretation 
to (3-7). The only terms in (3-7) that will give a non-zero 
contribution are those which contain the o,8 to an even 
power as the sum over states will cancel all terms which 
contain a o; to an odd power. With each (cj04) we associate 
a nearest neighbour bond (an edge of a graph) of the lattice, 
and since only even powers of the oj can occur in a given 
term, an even number of edges must meet at the ith vertex, 
or in the terminology of graph theory each vertex is of even 


degree. Thus with each term of (3-7) we can associate one 
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Ox more graphs, whose vertices are all of even degree. Such 
graphs, connected and separated, are referred to as non- 


magnetic graphs. 


The partition function can now be written as 
QN/2 9 qN/2 2 
Zy (8,0) = 2° (cosh K) {1+ 2 p(2,N) tanh K}, (3-8) 
L=1 

where p(2,N) denotes the total number of ways of construct- 
ing (embedding) in a given lattice all graphs of & lines 
whose vertices are all of even degree. The p(&,N) are in 
general polynomials of degree £ in N, hence the notation. 


N 


2 in (3-8) arises when we perform the sum over states. 


The procedure outlined in the preceding paragraphs 
is perhaps best illustrated by means of an example. The 
first six p(2,N) for the hypertriangular lattice are shown 


in Table 3.1 


In Table 3.1 we have enclosed in parenthesis 
those graphs contributing to a given p(%,N). This notation 
denotes the number of such graphs per lattice site and is 
commonly referred to as the lattice constant for a given 


graph. 


Instead of (3-8) what we are really interested in 


is the dimensionless Helmholtz free energy per site defined 
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TABLE 3.1 
The first six coefficients in the power series expansion 


for the zero field partition function of the hypertriangular 


lattice, 


p(l,N) = 0 p(4,N) = 0 
p(2,N) = 0 piS.N) =n (()) 
p(3,N) =N (A) = 3N 
= p(6,N) =N (4) ) +N (Pq ) 
+n(AA ) 
= 1/2 N2+141/2N 
TABLE 3.2 


The first four coefficients, calculated by the direct 


method, in the susceptibility series for the hypertriangular 


lattice. 


b,) = (4) = 3 b3 = (/V ) = 72 
be = (A) = 15 by = (AA) + (D- ) 
+ (A/) 


= 345 
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by 
-£/kT = lim {ln Zy (8,0) }/N. (3-9) 
Neo 
Hence 
-~£/kT = ln 2 + (q/2) In (cosh K) 
qN/2 2 
+ lim {ln [1 + = p(2,N)v J}/N, (3-10) 
N+ l=] 


where v = tanh K. It is customary to denote the logarithm 
of the configurational partition function per site by Ind, 


2. 76. 


qn/2 Q 
InA.= {ln [1.4.2 p(£,N)v"l]I}AN, (3-11) 
Q=1 


hence 


-£/kT = ln 2 + (q/2) In (cosh K) + lim lnA. (3-12) 


N+ 
Upon formally working out the logarithm in (3-11) one finds 
din A= & p Vis (3-13) 

where ep is just the coefficient of N in p(2,N). 


For example, for the hypertriangular lattice we 


have (c. f. Table 3.1) 
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ia A = é In [1+p(3,N)v2+p(5,N)v5#p(6,N)vet-eo], (3-14) 


Expressing (AA\) in terms of connected lattice constants 
(see Appendix C, page 162) and expanding the logarithm on the 
right-hand side of (3-19) we have. 
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The coefficients in (3-15) are alllinear in configurational 
lattice constants and no configurations with cutting points 
enter, a result wnich exactly parallels the Mayer theory of 
cluster integral expansions in the theory of imperfect gases. 

A rigorous derivation of the cluster expansion method as 
applied to lattice statistical problems has been given by Sykes 


Sc al. (1966). 


The reason we must stipulate N*- is that if we do 
not, the partition function will just be a finite polynomial, 
which is everywhere analytic and hence would show no singular 
behavicur. Hence the free energy per siet for a regular 
lattice of coordination number q and N sites in the limit 
Neo can be written as 


oo 1 ; 
(=f/eT) = in oie “~e/2). ln (cooshn<K) + Py dh (3-16) 
R=1 ; 
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where we now have 


inher oe ye 


(317) 
Q=1 % 


Using the techniques outlined above we have 
derived the series expansion for the free energy of the 
hydrogen peroxide lattice up to order twenty-eight and 


for the hypertriangular lattice up to order fourteen. 


The dimensionless free energy per site for the 


hydrogen peroxide lattice is given by 
(p£/kT).= 1n.2 +, (3/2). 1n, (cosh K) +,.1.1/2 v'° 
Poe LO A/a ve + 23172 ve" 
+ 6 3/4 v?® + 178 1/2 v2? + 484 v2" 
+ 774 v7®& + 3023 1/4 v*® + eee, (3-18) 
and for the gece eseet ier lattice by 
(-£/kT) = In 2 + 3 ln (cosh K) + v® + 3v° + 14 1/2 v° 
+ 33v’? + 57v® + 218 1/3 oe el L/2" 07" 
+ 3045v?! + 8925 1/4 v'? + 32676v'® 


129352 1/2 wi* + eee , (3-19) 
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The coefficients of v in (3-19) up to and includ- 
ing v'* were calculated by graphical enumeration on the 
hypertriangular lattice. The last two terms of (3-19) 
were obtained from (3-18) via the star-triangle transfor- 


IMation as outlined in Chapter V. 


The configurational data required to derive the 
first twelve terms of the zero field partition function and 
dimensionless free energy for the hypertriangular lattice 
are given in Appendix A. The only A eraeved graphs which 
contribute to the first twenty-eight terms of the zero 
field partition function and dimensionless free energy for 
the hydrogen peroxide lattice are the p-graphs (polygons) 
listed in Appendix D, where we have tabulated separately 
for both lattices the "Star" graphs (see Sykes, Essam, 
Heap, and Hiley 1966), which are required in the calcula- 
tion of the lattice constants of the separated configura- 


tions listed in Appendix C. 


The remainder of the thermodynamic functions of 
interest, with the exception of the initial susceptibility, 
can be obtained from (3-16) by differentiation. The 


internal energy per site is given by 


9 
U = kT? a (-£/kT) (3=20) 
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and the zero field specific heat per site by 
C= 93U/oTt (3-21) 
For the hydrogen peroxide lattice we have 


U/J = - 11/2 v - 15v? + 15v!! - 2iv!3 - 147y!5 
- 254v!7 + 288v!9 - 3792v2! - 7689v23 


~ 8508v2> + 64527v27 + ooo, (3-22) 


C/k = K*2{1 1/2 - 1 1/2 v2 + 135v8 - 300v?9 + 438y!2 
+ 3932vi4 + 2130v!5 - 9807v!8 + 85104v29 
+ 97215v22 + 35853v24 + 1529529v2& +e0e, (3-23) 


whereas for the hypertriangular lattice we have 


U/J = - 3v = 3v2 - lav* - B87v5 - 216v® — 369v’ 
- 1734v® - 9159v9 - 31530v!9 - 97488v}! 
- 391293v!2 = '1703832v!3 - -e-, (3-24) 


C/k = K2[3 + 6v - 3v2 + 42v3 + 1248v° + 2148v® 
+ 12576v’ + 79848v®8 + 301428v2? + 989937v!° 
+ 4380216v!! + 21077448v!}2 + cee, (3-25) 


In terms of the variable K alone we have for the 


hydrogen peroxide lattice 
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C/k = 1.5K? - 1,.5K* + K® - 0.566666K®+ 135.2952381K!0 
= 660.1462434K!2 + 2002.070046K!4 - 2394.318471K!6 
+ 401,.5388552K!8 —- 8224.447893K2° + 145141.528K22 
= 677416.1035K2* + 1829923.386K2® - 1989842.638K28 


ofe yd (3-26) 
and for the hypertrinagular lattice we have 


C/k = 3K2 + 6K3 - 3K* + 40K5 + 437K® + 1206. 8K? 
+ 1566.866667K® + 10526.47619K? + 76074.59048K!° 
+ 274283.2423K!! + 781916.7583K}2 
+ 3515116.625K!3 + 18106355.75Ki+ + eee , (3-27) 


The entropy per site is obtained from the 


relation 


S/k = (U/kT)-(£/kT). (3-28) 
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B. The Initial Susceptibility 


In the previous section we have concerned our- 
selves with the calculation of the zero field partition 
function and the thermodynamic variables derivable from 
it. In order to derive the initial susceptibility per 


Site defined by 
Xo = CO" /0H") (£/kT) [a » (3-29) 


we must also consider the effect of the factor containing 


the magnetic field H in (3-3). The end result is given by 


X = xX,/(Nm?/kT) = 1+ 2 2 by vv, (3-30) 


where x is the reduced initial susceptibility per site. 

bp is the coefficient of N in the total number of ways of 
embedding in the lattice all graphs of £ lines and exactly 
two vertices of odd degree. Such graphs are referred to 

as Magnetic graphs. This method of calculating the initial 


susceptibility we shall refer to as the direct mdthod. 


The calculation of the first nine terms in the 
susceptibility series for the hypertriangular lattice using 
the Gerace method is given in Appendix E. For example, 
those graphs contributing to the first four terms of xX for 


the hypertriangular lattice are shown in Table 3.2. 
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The calculation of the initial susceptibility by 
the method of (3-30) requires the enumeration of graphs 
having vertices of degree one. As these graphs are 
exceedingly numerous , this method is very susceptible to 


error and excessively labourious. 


Sykes (1961) discovered a theorem which greatly 
reduces the labour involved in calculating the series for 
the initial susceptibility. Sykes' theorem, as later 
derived from graph-combinatorial arguments by Nagle and 


Temperley (1968) can be written in the form 


ae nee 2 2 => Q 
X = 1+ (l-ov) [qv(l-ov)-2(l-v*)= g(p) 2£v 


Pp 
z) > Q gq 
+ (ltv)? = g(p) v’ = s(s=2) py 
Ni s=0 
p 
+ 2 (1l+v) ? 5 gb) v* (s -1) (ss-1)] (3-31) 
> i J : 
Pp 


where © = (g-1), q being the coordination number of the 
lattice. p= (p teheDytyatPs 7 BEY P,) 48 a (qtl) tuple or 
vector partially describing a graph by the number P, 

of its vertices of degree s. g(p) is the number of graphs 
with description D, which may be embedded in the lattice 


and for which po = 0. 
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In (3-31) the singly primed summation is over 
all graphs (nonmagnetic) whose vertices are all of even 
degree, while the doubly primed summation is over all graphs 
(magnetic) with exactly two odd vertices, i and j, of odd 
degree S; and ss. Note that if either s; or Sj or both 


J 
equal one, the graph does not contribute to (3-31). 


Using Sykes' theorem we have derived the first 
twenty-seven terms in the expansion of the reduced initial 


susceptibility per site for the hydrogen peroxide lattice 


X = 1+ 3v + 6v* + lav* + 24v* + 48v° + 96v° 
+ 192v’ + 384v® + 768v? + 1506v'® + 2952v?? 
+ 5814v'? + 11448v'> + 22494v'* + 44192v'® 
+ 86514v'* + 169320v'’ + 330858v’* + 646524v"? 
+ 1263696v7° + 2469456v** + 4817616v"* 
+ 9397740v?> + 18315734v2" + 35694200v** 
+ 69538944v7® + 135463148v*? + oes (3-32) 


For the hypertriangular lattice we have obtained 
the first twelve terms in the series for the susceptibility 
by way of Sykes’ theorem. We have also added one more 
term by star-triangle transformation from (3-32) (c.f. eq'n 
(5=37)). The initial susceptibility for the hypertriangular 


lattice is given by 
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xX = 1+ 6v + 30v* + 144v* + 690v" + 3276v° 
+ 15360v® + 71532v’ + 332166v® + 1537836v° 
+ 7097244v*° + 32677584v** + 150203012v'? 
+ 689366094v'*® + oo» (3-33) 


The graphs required to derive the suscepti- 
bility series for the hypertriangular lattice via 
Sykes’ theorem are given in Appendix B. For the hydrogen 
peroxide lattice the only connected graphs required in 
calculating the first twenty-seven terms of the suscepti- 
bility in this way are the p-graphs and theta graphs of 


Appendix D, and the dumbell graphs of Appendix J. 
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CHAPTER IV 


LOW TEMPERATURE SERIES. EXPANSIONS OF 


FERROMAGNETS AND ANTIFERROMAGNETS 


A. The Direct. Method. 


Consider an Ising model. of a ferromagnet at 
T = 0- All the spins are aligned and give rise to a 
spontaneous magnetization. As the temperature is raised 
above T = 0 thermal fluctuations will perturb the. ground 
state. The probability of any given perturbation will be 
given by the appropriate Boltzmann factor. As overturning 
almost any spin causes an increase in energy, the most. 
important perturbations will correspond to a relatively 
few overturned spins. Thus we. can group the perturbations 
according to the number. of overturned spins, the energy 
of any particular perturbation depending on the relative 


positions of these spins. 


Denote by -J the interaction energy of a pair 
of aligned first neighbour spins and by -mH the interac- 
tion energy of a spin aligned with the magnetic field H. 
Then if we overturn s spins of magnetic moment m with 
r first neighbour bonds between them there will be a 


gain in energy of 
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2(qs-2r)J3 + 2msH , (4-1) 


where qi? sis. the coordination number of the lattice, i.e. : 
the number of first nearest neighbours. Thus the 
Boltzmann factor for a perturbation having its energy 


given by (4-1) will be 
exp {[-2(qs-2r)J-2msH]/kT}. . (4-2) 


Letting exp (-43/kT) "= u.= 2% and, exp(=2mH/kT)-= 1. , 


(4-2) can be written as 


1 
2d Sar ue 


exp{[-2(qs-2r)J-2msH]/kT} = u (4-3) 
The energy of the ground state is given by 
e. = - n(dqo+mH) (4-4) 
fe) 24 ‘i 
Thus we can write the partition function as follows 
Qy(u,H) = exp(-e,/kT) Ay (u,u) + (4=9) 
where 
Zqs-x Ss 
Ay (uu) =99) Gy (err) Th a (4-6) 
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For a given lattice the G.(s,r) in (4-6) denote 


n 
the total number. of strong embeddings (c.f. Domb 1960, Sykes 
et al. 1966) of s spins with r . nearest neighbour bonds 
between them. In general Gy ls.X). will be a polynomial in 


N. For the hydrogen peroxide. lattice the first few 


Gy fS rE) are given by 


Gian ple = 14N 
vw 2 
G,,(2,0) = ty 2-2Nn 
Neen Z, 
Gy (3,2) = 3N (4-7) 
af (syih = 1iN?-9N 
1 ete! 2 
i 1 Sh 2 i 
G, (3,0) =_ gN 2N +65N 
The free energy per spin is given by 
eet Satin im. SE ag SO) (4-8) 
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Exactly the same arguments are valid here as were applied 
in Sect. A of Chapter III in dealing with the taking of 
the limit, i.e. upon taking ar [lin Ay (au) J/N we find 

=> OO 


that only the coefficient of N in Gy(s,r) fails to vanish 


and hence 


. *QS-r S$ 
lim Cin A (usu) I/N z# ln AC(uyy) = ¢£ C[s,rJu 


N-ve s,r (4-10) 
where by [s,r] we mean coefficient of N in Gy (sr). 


It is customary to group the expansion (4-10) 
as a series in powers of yp, the successive coefficients 


then being finite polynomials in u, L(u), defined by 
In A (u,n) © © Le(u)u’. (4-11) 
s 


The polynomials L,(u) are called low temperature ferro- 
magnetic polynomials. A given polynomial L, (u) is 
Bist ietely specified by a knowledge of all perturbations 
with s overturned spins, 

s(s-1) 


; 3qS—-P 
Beka) =r s [s,r]us?” : (4-12) 
r=Q 
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When using the method outlined above, which we. 
call the direct method, one. usually groups the perturbations 
according to the topology of their nearest neighbour 
linkages. As an example, we have listed in Fig. 4.1 those 
perturbations on the hydrogen peroxide lattice contributing 


to [7,6]. 


The low temperature ferromagnetic polynomials of 
the hypertriangular lattice which have been derived via 


the direct method are listed in Appendix PF. 


We will now describe the low temperature anti- 
ferromagnetic counting problem.. We consider a lattice which 
can be decomposed into two equivalent sublattices. By this 
we mean that each site of the lattice can be coloured in 
one of two ways such that a given site is surrounded by 
nearest neighbours all of the opposite colour. In the 
terminology of graph theory the lattice is bi-colourable 
(c.f. Fig. 2.6). Let one. of the sublattices be denoted by 
A and the other by B , and. to avoid the use of fractions 
we consider a lattice of 2N. sites, N of which are A-sites 
and N B-sites. Furthermore,. as. we are now discussing 


an antiferromagnet for which J < 0 , we shall take J” = -J. 


For small fields the lowest energy state is now 


one of antiparallel ordering in which each A-spin is in 
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Fig. 4.1. Topological breakdown of the 
perturbations on the hydrogen peroxide 
lattice contributing to [7,6]. The 


coefficient of z* in L,(z) = 143. 
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a state which is the opposite of its first nearest neighbours 


each of which is a B-spin.. If we perturb such a state by 


Overturning s A-spins. and t- B-spins with r nearest 
neighbour bonds between them the resultant gain in energy 


will be 


2[q(st+t)-2r]J°+2msH-2mtH , (4-13) 


where we have chosen the direction of the field as that of 
the unperturbed A-sites. It: should be pointed out. that 
equation (4-13) for the case of zero field reduces to 
(4-1), which simply expresses the well known result that in 
zero field the configurational problem for a loose packed 


Ising antiferromagnet is isomorphic with the ferromagnetic 


one. 
We now introduce the variables y = exp(-2J°/kT) 

and W = y*.. The energy of the ground state is: given: by 

a = -NqJ /2. Hence the partition: function is given by 
a a a 
Qy (wri) = exp(-e /kT) Ay (wu) , (4-14) 
i 
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The free energy per spin is given by 


iligaer | 
# belweqar Swen etnids he(wlanosed (4-16) 


(stt)-r s -t 
oe 


a 
where In A*(w,u) = : [stt,r]w2? (4-17) 


Sitar 


For zero field the low temperature expansion of 


the partition function will be given by 


J t)- 
In. A? (w;1)-= ) [s+t,r]w 2d (stt) hie (4-18) 
s,t,xr 


where as in the discussion of the ferromagnetic problem 
we denote by [s+t,r] the coefficient of N in 


Gy (sttix) 


Now consider the case of a ferromagnetic loose 
packed lattice. Again we denote one sublattice by A and 
the other by B. To the A-sites we assign a magnetic 
moment m and to the B-sites we assign a magnetic moment 


A 
mM, - In our previous description of the ferromagnetic 
problem we were concerned with perturbations of s_ spins 
with r first neighbour bonds between them. We are now 


concerned with enumerating perturbations consisting of 


s A-spins, t- B-spins , with r first. neighbour bonds 
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between them. The coefficient of N in the total number of 
ways of embedding: s. A-spins, +t. B-spins, with. r..first 


neighbour bonds between them will be denoted by ([s,t;r] 


The expression for the: free energy per site will 


be given by 


Le ac ee ae | | , . 
fo =. 599 mH m,,H REO dns ia a his (4-19) 


where 


us aed ace (4-20) 


3 
Inwhs(unbavie= fe, le xtsr)u27 


The variables wu and vy are defined by 


exp (-2m,H/kT) 


I 
i= 


exp (~2m,H/kT) 


I 
<E 


If we set v= 1/u and u=w we have 


a ES 5 
in. A (w,y) = ) ispeortue ag gat ne (4-21) 


Sop ti, Ic 


which can be looked upon as a formal isomorphism with (4-17). 


We have retained the [s,t;r] of the ferromagnetic problem 
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even in the presence of a field. The isomorphism between 
& 

Ree kh (4-17) cand. Linsey) (4=21)-) is: enly. formal 

since the two series have different regions of physical 


validity. 


The low temperature polynomials for the sub- 
lattice ferromagnetic problem are defined by 
(stt)-r 


a 
NE; =) [s,;tzr]u2? 


(4-22) 
Sit bs 


Again the method of topological breakdown can be used. to 


enumerate the various [s,t;r]. 


While the above.method (the direct method) is 
simple and quite straight forward it is difficult to derive 
more than a few terms as the number of perturbations 
contributing to a given .[s,r] or [s,t;r] increases 
rapidly. A further disadvantage is that there is no 
systematic way of checking whether one in fact has listed 
all the perturbations. correctly or. if one has calculated 
the contribution of a given perturbation correctly. 
However, there is another method, being both more 
sophisticated and more powerful and essentially self- 
checking, by which we. can attack the problem of enumerating 
the required perturbations. We shall proceed to describe 


this method in the following section. 
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B. The Shadow Lattice Method. 


We shall now describe a method, which was introduced 
by Sykes, Essam, and Gaunt (1965), of enumerating low 
temperature perturbations on a loose packed lattice. This 
method eliminates the need for any topological description 
@f-the “perturbations contributing to a*given “[s;t:r] and 

instead simply provides us with the numerical values of 


the required {[s,t;r] 


We define a generating function F by 


F(X,Y,b) = J. [s,tsr]xSy*p™ . (4-23) 


Ste 


If we could find a closed form expansion for F we would 
in fact have the equivalent of an exact solution to the 
low temperature enumeration problem. Instead of using 
the generating function .F., we shall find it more 
convenient to work with a partial generating function 


FB. (X,b) defined by 


F(X,Y¥,b) = } v\F, (X,b) , (4-24) 
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F; (Xibk-se J gisvdarlx bp. (4-25) 
S75 


The partial generating function F, (X,b) is equivalent to 
an exact solution of the problem when the number of over- 
turned spins, A., on one sublattice is held constant. In 
fact Since the sublattices are symmetric, a knowledge of 
the first n partial generating functions makes it 


pessible. to determine’ all the [s,t;r] for s.+ ¢ < 2n +71. 


We shall now show how one can derive the first 
three paratial generating functions for the hydrogen 
peroxide lattice. The derivation of partial generating 


functions of higher order will be quite straight forward. 


The zeroth. order. generating function Fo 
corresponds to configurations for which all the sites are 
on one sublattice, say. the A-sublattice . As the sites of 
a given sublattice are not nearest neighbours of each 
other, r = 0 always. Since we have an infinite lattice of 
2N sites, the first site on the A_ sublattice may be. 


chosen in N ways, the second in (N-1) ways and so on, 


If we want to choose s A-sites andno  B-sites 


the number of ways we can do this is simply 


N(N-1) (N-2)-**(N-s+l)/s° . (4-26) 
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We want’ [s,0;0], which. is simply the coefficient of N. in 


(4-26), hence 


[s,0;0] = (-1)8t+ (4-27) 
s 
Therefore for Bo we have 
F(X ,b) es In(i+x) . (4-28) 


To derive Be we must choose one site on the 
B sublattice, and there are N ways of doing this. Upon. 
choosing the B-site we observe that it casts a "shadow" 
on its three nearest neighbours, all of which are A-sites 
(see Fig. 2.6 and Fig. 4.2(a)) in the sense that if an 
A-site is chosen from. the sites of the shadow a nearest 


neighbour bond is formed. 


If we choose a A-sites from the three sites. 
in the shadow and §8 . A-sites from the remaining (N-3) 
Sites we shall obtain a nearest neighbour bonds in 


N (3) (N-3) (N-4) +++ (N-B+1) /B! (4-29) 


ways. The coefficient.of N in (4-29) is given by 
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Fig. 4.2 (a) A Besite (0) and its corres- 
ponding "shadow" composed of three A=-sites 
(e). The dashed lines indicate nearest 
neighbour bonds while the solid lines con- 
necting the A-sites form the "shadow." 

(b) The possible overlappings of the 
shadow cast by two B-sites. The figures 
refer to the hydrogen peroxide lattice. 

The A=-sites (and also the Basiites) each 
form the sites of a hypertriangular lattice. 
The occurrence factors are just the number 
of strong embeddings of the given configura- 


tion on the hypertriangular lattice. 
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Shadow Occurrence Factor 
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B 
v (3) ; (4-30) 


The term in ae corresponding to (4-30) is thus 


=i 


g) ’ (4-31) 


Cr x a 
fo} 
and therefore by the binomial theorem we have 


F (X,b) =. (1+bX)2(14x)7 °°. (4-32) 


To derive. F .we turn over two. B=-sites and 
2 


observe that the shadows may overlap as shown in Fig. 4.2(b). 


Consider. the configuration of two triangles 
placed corner to corner. in Fig. 4.2(b). Since the two 
triangles have one corner. in common, only five A-sites 
are shaded by the two B-sites . Of these five sites the 
one on the common. corner will create two first neighbour 
bonds if chosen since it is adjacent to both B-sites 
If one of the remaining four A-sites is chosen-one first 
neighbour bond will be formed. If an A-site other than 


one of these five is chosen no bonds will be formed. 


Let us now. choose from this configuration a 


A-sites , which each form two first neighbour bonds, and 
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8B A-sites , which each form one first neighbour bond, and 
yY A-sites' which form.no first neighbour bonds. The 


ao A-sites can be chosen. in 


ways; the 8 A-sites can be chosen in 


DBD & 


ways, and the remaining y A-sites can be chosen in 
(N-5) (N-6) °°* (N-5-y+l)/y! (4-33) 


ways. 


The corner to corner configuration of two 
triangles can be placed on the hypertriangular lattice in 
3N ways, therefore on the hydrogen peroxide lattice the 
number of ways of choosing two B-sites casting a shadow: 
of the corner to corner type plus a+f+y A-sites with 


2a + 8 bonds between them is 


3n (2) (2) (N=5) (N-6) *** (N-5-y+1)/7! (4-34) 
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and the term linear in’ N+ ‘in: \this-.is simply 
ie A eS 
3 (4) (a) su ; (4-35) 


Hence the centribution..of configuration. (1) of Fig. 4.2(b) 
to F is given by 
a 


ae cee 
a,B8,Y . 


paDexi Gy (exe yx! = 3(1+b2Xx) (1+bx) * 


oe CLE 


(4-36) 


For the two separated triangles there are six 
A-sites that can form a first neighbour bond with a 
B-site and hence (N-6) A-sites that will form no first 
neighbour bonds. Choosing a A-sites. from the six 
A-sites of the shadow and Bg A-sites from the remaining 
(N-6) A-sites, we can form on the hydrogen peroxide 
lattice configurations. of a+ 8 A-sites and two B-sites 


with a first neighbour bonds in 


(gn?-35N) (8) (N-6) (N-7) *** (N-6-8+1) /8! (4-37) 


ways. The: coefficient of N in. (4-37) is given by 


ipa ih i ae 
rae he be 


‘ ie 
COT ae ek oro 


‘(edety) (Ks deloe = Yea i yy ety GY” cna cf) i 
ya 


2 (ae~s) a P oe Aid A - j 


; ; | te 


xie Sue evens eelensiad beseisuse ad edd so% , in 


! n 
7 ry! 7" i 


& tHtiw nod wwodripter sad “& MOT ee iy ott 
sesnt on mot Litw reds aothe-4 alles ant | 


—. 


‘ iii pit moxt serterh # paiscoda 
gion edmbi ett mad aad aA! 8 fees \wabadte 9 4 ba- 
4 meee: mi ‘ ae ; 


| aliesteorseg insigitlge oy me MCR" £18 ibe a ni es 


is rh ue 


| dialed ond bas” seate-n ae + es reiuphiaes soLsI6 
ih ; ; nie v i ay py 

i } | my: aa ag ‘thot eee BYP % Bab Le 

| | . ie ot pe arth ar i 


{ i ; ¥ i Nie ; 7 ’ | : r \ VA : 
; i ‘ ve i bul : Bag S re ' ole ore vP wy ea 


es eee pane tn age ty 


i . f | | i : Fi: : mn | M a aA aie ‘ i } i" we) AA 
a | I ky ‘DE bane bee ail THM ; mn ee 4 1) 
Aiwa) saab 
Gant 
i ae ane Fost a ii an 


83 


3h (6 (78 


pop eR) 
Hence the separated. configuration gives the following 
contribution to F 
= 1 6 =6 
= - 35 (1+bx) (1+X) 

(4-38) 


Therefore the second partial generating function for the 


hydrogen peroxide lattice is simply 


F (X,b) = 3(1+bx) * (1+b*x) (1+x) 5 


ak 


- 35 (1+¢bx) (14x) 7 : (4-39) 


The general term in FL (X,b) will be of the 


form 


atB+ytee° (4-40) 


(1+bx) % (14b2x) 8 (14b 3x) Ve + / (14x) 


which we shall denote by the following code 


UNG Oy Be Via ty A= OF Bey ee (4-41) 
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which forms a convenient déscription for supplying such 
data to a computer. The brackets are closed when the 
factors to the right terminate. However, any intermediate 


zeros must be included in the code. 


‘Very simply the code tells us the following. 
For a given "shadow" composed of d A-sites, a of these 
each form one first neighbour bond with the B-sites, 
8B A-sites each form two first neighbour bonds with the 
B-sites, and so on. The coded form of F, and F, are 


shown below. 


F,(X,b) = (3,3) 
° (4-42) 


8 
F, (X,b) = 3 (S74 ASS ters) 


In this manner we have calculated completely 
the first seven partial generating functions for the 
hydrogen peroxide lattice and we have also calculated 
the leading terms Of Fag, For Fyp, and Fi, which are 
listed in coded form in Appendix G. The shadow con- 
figurations required to calculate the partial generating 
functions of the hydrogen peroxide lattice are listed 


in Appendix H. 
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In order. to.obtain the various [s,\;r] from the 
partial: generating functions it is necessary to expand 
them in the form of. a power series. This can only 
conveniently be done:..on a computer and for this purpose. 
the coded partial generating functions form the input data. 
We have written a program to perform such a task. However, 
as the output is. very extensive we shall not list it here. 
We have pointed out: previously that this method is 
essentially self-checking since the two sublattices are 
symmetric, ise@éidserA; ud =adAeserd., This  isea»great 


advantage of the shadow lattice method over others. 


As an example of how we use the partial generating 
functions consider the topological description of the 
perturbation of Fig. (4.1). We need the coefficient of 
b® for all s+t = 7. Expanding the first three partial 
generating functions we find the coefficient of b® is 


given by 
SXVe peausOxe Vil) & AVAOKEY ot aS? 1. (4-43) 
Adding the coefficients together we obtain 
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which is the same as [7,6] quoted in Fig. 4.1. 


The first seven partial generating functions are 
sufficient to determine completely the first fifteen low 
temperature ferromagnetic polynomials for the hydrogen 
peroxide lattices: We have siss calculated the leading terms 
of higher order polynomials up to and including Loo (2). 

The low temperature ferromagnetic polynomials for the 


hydrogen peroxide lattice are listed in Appendix I. 


In addition to the energy, specific heat, and 
initial susceptibility, which have been previously defined, 
we are also interested in the spontaneous magnetization 


defined by 


ie Sa 
3H Heo. (4-44) 


From the polynomials listed in Appendix I, we have calcu- 
lated ferromagnetic low temperature series expansions for 
the Ising model on the hydrogen peroxide lattice. The 


series are as follows: 
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(-£/kT) = (q/2) K + 23 + 11/2 24 + 3 25 + 5 26 


+9 27 +15 3/4 28 +128 1/3 29 + 52 1/2 210 


"+ 108 z!! + 254 212 4+ 657 2}3 


+°1738-1/2° 22" + 4511 1/5275 


+ 11463 3/8 zi& + coe, (4-45) 


configurational energy 


U/I = (=q/2) + 6 23 + 12 z¥ + 30 25 + 60 26 


oo 
- 
” 


+ 
specific heat 


C/k(lnz)? = 9 


+ 


+ 


op 


1426 27 + 252 25:+ 510 29 + 1050 2} 


2376 z!} + 6096 z!2 + 17082 2}3 
48678 zi* + 135336 z!5 + 366828 z!6 
nee (4-46) 
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z3> + 24 z4 + 75 25 + 180 26 + 430 1/2 2? 

1008 28 + 2295 29 + 5250 z!9 + 13068 z?!! 

36576 ge 4 111033 2/3 + 347462 21% 

1015020 2!5 + 2934624 z!& + oo, (4-47) 
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spontaneous magnetization . 


I =l=2 23 = 6 z+ = 18 25 = 48 2& - 126 27 — 234 28 

- 830 29 = 2154 z!9 - 5784 z!! - 16146 z}2 

- 46302 z!3 - 134082 2)" — 387724 215 - 1119024 zJ6 

pod. te, (4-48) 


and initial susceptibility 


¥ = 4 23 + 24 z* + 108 z5 + 416 z& + 1476-27 + 4968 28 
+ 17156 z9 + 51912 2)° + 167568 2}! + 545456 z}2 
& 1787364 zi3 + 5838360. 2!* +.18949712 2!° 
+ 61171104 zi& + eos , (4-49)° 


The ferromagnetic low temperature series expansions 


for the hypertriangular lattice are listed in Appendix Fe 
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CHAPTER V 


ISING MODEL TRANSFORMATIONS 


A. Zhe Star-Triangle Transformation 


The star-triangle transformation was originally 
developed by Onsager Le Wannier 1945). By utilizing 
both the star-triangle and dual transformation Onsager 
was able to calculate exactly the critical points of the 
triangular and honeycomb lattices. In this work the dual 
transformation is of no direct interest as it is applicable 


to two dimensional lattices only. 


The star-triangle transformation replaces a star 

consisting of three spins interacting with a central spin 
by a triangle of spins interacting with one another 
(Fig. 5.1). The partition function for the star (Fig. 5.1(a)) 
is given by 

fy = & exp Ky (ago ,/4 o995 + 699) (5-1) 

o9o=tl 

where Ky = Jy/kT. Since the o variables commute and obey 


the Van der Waerden identity (3-5), we.can write (5-1) as 
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Zy = cosh * io Pir ie + (050, + G90, +0,90,) tanh Ky 
yan 


+ “e 2 
(0,0, oo, + o,0,) tanh* K 


+ 9,0,0,0, tanh® K,]. (5-2) 
Upon performing the sum over %9 in (5-2) we obtain 


Z, = 2 cosh® ky Z [1 + (6,0, +0,0,+0,) tanh? K,] (5-3) 
O1,02,03=+1 


for the "star" partition function. 


The partition function for the triangle of spins 
(Fig. 5.1(b)) is given by 
01,02,03=+1 
where K, = J, /kT. Again we make use of the properties of 
the o variables, commutivity plus the Van der Waerden 
identity to obtain 


2, 5 cosh® K, £ [1 + (o:02 + 0103 + 0203) 
O,,02,0,=t1 


x (tanh K, + tanh? K,) + tanh? Kal, (5-5) 
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which can be cast into the form 
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Comparing (5-3) and (5-6) we see that if we 


choose K, to satisfy 


cosh K, sinh K, (cosh Ky + sinh Ky) A 
SS ee ee ee tanh ky 


weg ° 3 
ste K, + sinh* Ka (5-7) 


the partition functions of the star and triangle will be 


related to each other by the simple factor 


¢ = 2 cosh® K,/(cosh® K, + sinh® Ky), (5-8) 


A 


that is 


Equations (5-7) and (5-9) constitute the star-triangle 


transformation. 


The temperature variables (Ky = Jy/kT, Ky= J,/kT 
of star-triangle lattices are related through (5-7), which 
can be greatly simplified by the use of certain well known 


hyperbolic identities to yield 
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exp (4 K,) = 2 cosh (2 Ky) - 1. (5-10) 
In terms of the variables z = exp (-2J/kT), we have 

ZA = 2y/(l - zy + zy), (S-11) 
and in terms of v = tanh (J/kT) = tanh K, we have 

vy = (vy + vi)/(1 + vA). (5-12) 


If we apply the star-triangle transformation to a 
lattice of connected stars we can write 


N 
Z (Ky) = & Z 


- ee (5-13) 


N 


where Y and A denote "star" and "triangle" lattices respec- 


tively. The subscripts 2N and N refer to the number of sites 


on the respective lattices. In three dimensions we have 


y HP (hydrogen peroxide) 


A HT (hypertriangular) 


and in two dimensions we have 


K 
ii 


HC (honeycomb) 


> 
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T (triangular) 


The significance of the star-triangle transforma- 
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relationship between the critical points of the hydrogen 
peroxide and hypertriangular lattices, and hence provides 
us with a valuable check of our numerical analysis of the 
series of Chapter III. The star-triangle transformation 
also enables 42 to derive from the hydrogen peroxide 
lattice the partition function and intial susceptibility 
of the hypertriangular lattice. We shall have more to say 


about this in Section C of this Chapter. 
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B. The Decoration Transformation 


The decoration transformation was first suggested 
by Syozi (1951) and later developed by Naya (1954). The 
decoration transformation combined with the star-triangle 
transformation enables one to calculate exactly the 
solution of the Ising model on the kagomé lattice from the 
solution of the Ising model on the honeycomb lattice. In 
this work we shall combine the decoration transformation 
with the star-triangle transformation to calculate the pro- 
perties of the Ising model on the hyperkagomé lattice 
(Fig. 2.9) from the properties of the Ising model on the 


hydrogen peroxide lattice (Fig. 2.1). 


The decoration transformation replaces a central 
Spin interacting with two neighbouring spins plus the 
interaction with the magnetic field by a single bond con- 
necting the two outer spins plus a new interaction with the 
Magnetic field. The partition function for the system shown 
mon Fig. 5.2(a) fe given by 


Zz. = i exp (Kp) 590; + Kp 5952 + Lp Og), (5-14) 
D Oo=+1 


oO 1 pO 2=+1 


where Ky = Jp/kT, Lp = mH,/KT . The subscript D denotes a 


decorated system. Again using the commutivity of the o 
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variables plus the Van der Waerden identity, (5-14) can 
be written as 
Z. = cosh? Ky) cosh lp £ {1 + (0; + og). tanh K 


y O_=tl 
0,90 ,=t1 


D tanh Ly 


+ 03;02 tanh? Kp + do [(o, + 92) tanh Kp + tanh L, 


+ 9,0, tanh® K_ tanh Lp]}. (5=15) 


a D 


Performing the sum over o,.we are left with 


Z. = 2 cosh? K 


D cosh Lp = [1 + (9, + o,) tanh Kp tanh Lp 


D 
O, pF gatl 


+ 6.0, tanh? KI: (5-16) 


l~ 2 


The partition function for the system shown in 


Fig. 5.2(b) is given by 


Z= ft exp (K 0,0, +Lo, +La,), (5-17) 


where K = J/kT and L = mH/kT. Using arguments previously 


given (5-17) can be written in the form 


sale 
. | cosh K cosh = | rales. ten es) 
1 + tanh K tanh* L Oy 02541 


, {tanh L_+ tanh K tanh L) , , |. |tanh K + tanh? L 
1 + tanh K tanh? L ? l1 + tanh K tanh? Li’° 


(5-18) 
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Comparing (5-16) with (5-18) we see that if we require 


tanh K, tanh L, = tanh L + tanh K tanh L 
1 + tanh K tanh? L 


(3-19) 


tanh” Ky = tanh K + tanh* L 
1 + tanh K tanh? L 


the partition functions for the decorated system (Fig. 5.2(a)) 
and the undecorated system (Fig. 5.2(b)) will be related in 


the following way: 


6. = 62, (5=20) 
where 


cosh? Kp cosh Ly 
Es 2 |———————-/|_ (1 + tanh K tanh’ L). (5-21) 


cosh K cosh? L 


From (5-19) we see that for zero field the rela- 
tionship between the temperature variables (Kp = Jp/kT, 
K = J/kT) of decorated and undecorated Ising lattices is 
given by 
tanh* Kp = tanh K (5-22) 
or 
5 = V. (5-23) 


In terms of the variable Zz we have 
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z= 2 Zp/(1 + 2°). (5-24) 


If we apply the decoration transformation 
(5-20,-21) to a lattice of W sites and coordination number 
q we will get a decorated lattice of W+ (gW/2) sites. The 
corresponding partition functions will be related by 

2 (Kp) tle inte (5-25) 

N+ (Qv/2) N 

Applying the decoration transformation to the 
hydrogen peroxide lattice (Fig. 2.1) we get the decorated 
hydrogen peroxide lattice (Fig. 2.8). The partition func- 
tions of hydrogen peroxide and decorated hydrogen peroxide 


lattices are related by 


3N ; 
Z (Kpyp) = & 2 (Kyp), (5-26) 
a HP Bae 


where DHP and HP denote the decorated hydrogen peroxide and 
hydrogen peroxide lattices respectively. We have taken the 
number of sites on the hydrogen peroxide lattice as 2N, and 
since the hydrogen peroxide lattice has a coordination 
number of three, we obtain 5N sites for the decorated 


hydrogen peroxide lattice. 


If we now apply the star-triangle transformation 
to the open circles (0) of the decorated hydrogen peroxide 


lattice (Fig. 2.8) we get the hyperkagomé lattice of Fig. 2.9. 
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The partition functions. for these two lattices. are related 


by 


Z (KK...) = c2Ng (Kx) 5=27 
SM, BHP Sn HK?’ ( ) 


where HK denotes the hyperkagomé lattice. The hyperkagomé 
lattice has only 3N sites as 2N sites of the decorated 
hydrogen peroxide lattice are removed in the process of 


carrying out the star-triangle transformation. 


The relationship between the partition functions 
of the hyperkagomé and hydrogen peroxide lattices is 
obtained by combining (3-26) and (5-27). For zero magnetic 


field we have 


N 
3 
Z (Kuz) =| 2 BRAN (5-28) 
3N_ eosh™ Eup 2N 


Equation (5-28) enables us to calculate from the hydrogen 
peroxide lattice the zero field partition function, engery, 


and specific heat of the hyperkagomé lattice. 


The relationship between the temperature variables 


(Ku = Jyx/kT, Kyp = Jyp/kT) of hyperkagomé and hydrogen 


HK 
peroxide lattices is obtained from (5=7) and (5-22). The 


result is 


exp (4Kpy) = (1 + 3 tanh Kyp)/(1 + tanh Kp). (5-29) 
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In terms of the variable v we have 


on 2 3 
ve (Vix + Var) / 4+ Vi) ¢ (5=30) 


and in terms of z we have 


2 


25K = Zup/ (2 a Zp) ° (5-31) 


The development we have outlined above has been 
concerned only with the transformation from the hydrogen 
peroxide lattice to the hyperkagomé lattice, which are 
three dimensional lattices. Exactly the same procedure 
holds for the transformation from the honeycomb lattice to 
the kagomé lattice, which are two dimensional lattices. The 
only change that is required is one of notation, HP + HC, 
and HK + K, where HC and K stand for the honeycomb and 
kagomé lattices respectively.. A very thorough review of 
the star-triangle and decoration transformations with 


arbitrary interactions has been given by Fisher (1959(b)). 
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C, Exact Relationships Between the Thermodynamic 


Variables 


The star-triangle and decoration transformations 
discussed in the previous section make it possible to 
determine from the hydrogen peroxide lattice the ther- 
‘modynamic variables (partition function, energy, specific 


heat, etc.) of the hypertriangular and hyperkagomé lattices. 


The partition functions of star-triangle lattices 
are related by (5-13), and hence the relationship between 
the dimensionless Helmholtz free energy per site for star- 


triangle lattices is given by 
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where (c.f. eq'ns 3-9, -10) 
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(3-13). From (5-8) we have 


in ¢ = ln 2 + 3 In (cosh KL) - 3 ln (cosh K,) 


- ln (1 + ve (5-35) 
Hence 
(Lye 3 Ay CE 
ae Ma ins (lst 7 +2 Pa Vy» (5-36) 


‘where the relationship between Vy and v, is given by (5-12). 


A 


From (5-36) we see that if we know the first 2N 
terms of the free energy for the star lattice (denoted by 
Y), we can calculate the first N terms of the free energy 
for the triangle lattice (denoted by 4). The same holds 
true for the internal energy and specific heat as they are 


found from the free energy by differentiation. 


In this way we obtained from the series for the 
free energy of the hydrogen peroxide lattice, for which we 
have calculated the first twenty-eight terms, the first 
fourteen terms in the series for the free energy of the 
hypertriangular lattice. The first twelve terms in the 
wéties fbr the free energy of the hypertriangular lattice 

were also calculated directly from the hypertriangular 


lattice (see Appendix A) and as far.as they go the two 
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series agree exactly. 


Fisher (1959(b)) has shown how one can relate 
the initial susceptibilities of star-triangle lattices. 


The result is simply 
xX (vy) = 1/2 [x (+v,) ae Oe Gn ee (5-37) 


As was the case with the free energy, the first 2N terms 
in the series for the initial susceptibility for the star 
lattice give rise to the first N terms in the seta: 
the initial susceptibility of the triangle lattice. Thus 
from the first twenty-six terms in the susceptibility 
series for the hydrogen peroxide lattice we obtained the 
first thirteen terms in the susceptibility series for the 
hypertriangular lattice. The first twelve terms in the 
series for the initial susceptibility of the hypertriangular. 
lattice were also calculated independently of the hydrogen 
peroxide lattice (see Appendix B). The results of the two 
methods agree exactly, indicating that for both lattices 


the graphical enumeration is correct. 


As was pointed out in Chapter I, in the critical 
region the specific heat and initial susceptibility can be 


represented by the asymptotic forms 
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x oo, OF 
C= A, (1 - K/K,) . 5-80) 


Pent gh tas KUK. dey (5-39) 


where A, and CL denote the amplitude of the specific heat 
and initial susceptibility respectively. By the use of 
equation (5-13) and (5-28) we can derive exact expressions, 
valid in the critical region, which relate the amplitude 


on one lattice to that of the other lattice. 


For example, for star-triangle lattices the 
partition functions are related by (5-13). Since we are 
rested only in the singular part of the partition 
function the factor 7 in (5-13) can be neglected. Thus 
for the singular part of the dimensionless Helmholtz free 


energy per site we have 


A 
(-e/et)” = 4 in 24 (x4) = ofet an zt. (KY)] = 2(-£/ur)*. 
N N EN 2N \ 
| (5-40) 
The specific heat is defined by 
C = kK2 92 (=f/kT)/2K2 (5-41) 
From (5-10) we have 
ex (~4K*) pee (sinh ot) : (5-42) 
Riprinng/’ espa oe aKy 


Hence in the critical region we have 
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2 


(KG)? exp (-8Kq) C°(K) = 2(KS)2 (sinh 2Ke)* (KE) (5-43) 


Upon substituting in (5-43) the asymptotic form (5-38) 
we get 

y a 
(K,)? exp(-8Kg) af (2 - K‘7K4 )7° 


i Ay2 Ae Grey Rane 4 Y¥,-G 
= 2(K )*(sinh 2K.) A, (1 ~ K./K) (5-44) 


Taking the logarithm of (5-44) and applying (5-42) we 


find 
Ke (~4K4) ( A dy? KA (ginn 2Kt)F? ca 2 x ety? 
o exp (-4KQ)(1 - K"/Ke) = Kg. (sinh 2K,) = K /Ko) 
(5-45) 
From (5-44) and (5-45) we find for the final result 
Yo, a Lay eins Y.2-a 
A/a, = 5 [(KQ/K,) exp (-4K’) sinh 2K] (5-46) 


which is valid for star-triangle lattices of arbitary 


dimension. 


An expression similar to (5-46) can be derived 
for the amplitude of the Ne cepeinaety of star-triangle 
lattices. For example, on a star-lattice, which has a 
coordination number of three and 2N sites consisting of 


primed and unprimed sublattices, the reduced initial 
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Beusceptibility per site is fiven by 
_ Y Roe. 3 | 
= xo/(2Nm2/kT) = = 2 (<oyo,> + 2 <ayo!> 
i=l z “ 
<o10,>). (5-47) 
Mmaese three sums are in gernal distinct, but in the critical 


region where the spin correlations become long range, the 


sums become equal, hence 


N v 
alee ] 4 Y 
xX = commas = Woy bk )/Zon(K ) (5-48) 
where 
N 
yk) = 2 ee Tee TS uae a eee AE 
i=l o,=t1 o =+)] Gaal oy=t a 
ae is — 
ol San 
x exp K £ £ aiag). (5-49) 
{=z] 621 


The subs*ript 6 indicates the three nearest neighbours of 
the ith spin. The star-triangle relation can be applied to 


(5-49) and we end up with 
xY(K") = 2 x4(K4), (5-50) 


Now proceeding exactly as in the case for the specific 


heat we find 
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Y 
oe 


/c) = 20(Ko/KQ) exp (-4K9) sinh 2x%]7’, (5-521) 


A 
+ 
Results similar to (5-46) and (5-51) can be 

obtained for the hydrogen peroxide and hyperkagome lattices 
(honeycomb and kagome lattices in two dimensions). The 
procedure is exactl; as outlined for the star-triangle 
lattices except now one has the intermediate decoration 
process. Instead of (5-13) one now uses (5-28) and in 


place of (5-42) one has 


HK 

exp (-4K, ) ——— 
S 9 HK 

= [(1 ~ tanh KEP)/(1 + tanh Ke) ap up: (5-52) 


The specific heat amplitude relation is given 


by 

HK HK age 

HP 
: HK r 

A+ 2 2 |f Kee) exp (uke) += #8nh 4e. (5-53) 

HP ; P 1 + tanh K 

A 3 Ke e 

ofp 


while the amplitudes of the initial susceptibilities are 


related by 

~¥ 
HK HK | a ~54) 
Gyo 3 aa lheg | (ugHK) 1 _- tanh Ke (5-54) 
HP ¢ HP 1 + tanh x 
Cy Cc 
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For the two dimensional honeycomb and kagome lattices 
one has only to make the following change in notation in 


(5-53) and (5-54): HP+HC, HK+K. 


The symmetry between (5-46), (5-51) and (5-53), 
(5-54) is quite obvious. In fact, having derived the 
amplitude relationships for the star-triangle lattices 
one could write down by inspection the expressions (5-53) 
and (5-54). However, (5-53) and (5-54) were derived 


independently of the results for the star-triangle lattices. 


The validity of the expressions (5-46, -51, -53, 
-54)can be tested for two dimensions. In two dimensions, 
where the star lattice is the honeycomb lattice and the 


triangle lattice is the triangular lattice, a = 0, 


y = 7/4, exp (4KE) = 3, and sinh 2ku" = v3 yield 
is (Oe ee ge 
Be Are 0.057994 


and 


cHC oT 2 1,132234°°° 


qi 


in precise agreement with the ratios as quoted, for 


example, by Fisher (1967). For the kagome lattice 
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exp CAKE) = 3 + 273 and we find 


eG 
AL/A, = 1,004091°°° 


and 
ep/cHe = 1.048238>-- 


The first ratio agrees exactly with the ratio as 
computed from the individual specific heat amplitude 
estimates as quoted by Fisher (1967). However, using the 
estimates of the amplitudes of the susceptibility from 


Fisher's table we find 


KoRn , 
Ci /Cy = Ou. 72s ee: 


This 8% discrepency indicates either that our relation 
(5-54) is incorrect or that the estimate of the amplitude 
of the susceptibility of the kagome lattice is in error. 
The other three expressions for the amplitude ratios are 
clearly correct and the fourth, as already pointed out, 
could be written down by symmetry, and thus we believe 
(5-54) is correct. Sykes (private communication) has 
recently reexamined the susceptibility series for the 
kagome lattice and could find no error. It is the view of 


the author that this contoversy is still unresolved. 
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CHAPTER VI 


SELF-AVOIDING WALKS 


A. Relationship To Ising Model 


Self-avoiding walks are random walks with the 
restrictions that there can be no immediate reversals and 
no self intersections. Hence self-avoiding walks define 
a non-Markovian process and the statistical properties 
of such walks are not well known. The generating function 
for self-avoiding walks is define by 


C(x) = 1+ £ cx , (6-1) 


bs 


where c is the number of self-avoiding walks of length 


R 
£ on a lattice. The problem is to determine the asymptotic 


behaviour of the Cy F 


The relationship between the self-avoiding walk 
problem and the Ising model has been studied by 
Temperley (1956) and by Fisher and Sykes (1959). A more 
recent review has been given by Domb (1969). It is now 
generally recognized that the relationship between the two 
problems is very subtle and not entirely clear. Nonethe- 


less the self~avoiding walk problem provides a very good 


; € 
Wy 


, - ; oY 
7 re ee 
| | Se 
ei wrnetaccin hi Ral * ‘-e ‘ a, a 
Y oie at se) | 
LAW DERG qabvantat 5 
, < ‘i = Lit fy is 
- Bay | 


Leon pada _of qideno +#/ 


ay ; 


“olay, ova waldew pardiore bia 
5 Me 


swt of edt mae ore sont Bit, 


ea i as te) uy orp Slow aonen vanotipeaedas ® 


i a hips! 
stiragord jedise lease. ont bas nseooad ma tvo 
‘ 7 ite 2 


ae 


av .owemt tiew ton ots adden if 


eee 
an rae antted wt Ss iaianait re 
a : i i a 
| 1 aell 
[de ) ; "hy? 7 +h (4) \ dae 
‘ | Lwa 2 is 7 : 
4 a) : - ; 4a ne p ; i J 4 
| fel to easy sseiuasuntan, 20 eer et ‘eile ak i: 2 
a, Diy i 

. digncyae odd eaigheteh of af aeldorg edt 


: 
>] j .' ¥ 
2 ; ‘ 
. 
Uy -* , my * 
.- = ~ 7 J 
. . 7 ‘ 
st ig . Ari 

: ‘ t 


, j , diew path tove any et! ‘edd weaned aidanaiseter e6t 


es hers ng Pua ya boreal mead ‘pm Lebon. ame od 
| fe aa (oebty ‘epwye baa arian salary pees ek: 
Wh ; iiss | Wer vy ob ot , (@e@8) cose, et yapeseanty 
st at anes aang a Pad 
ace seta aslo ylexisae - 


| oll Fai. 2a ae 


— 


iy 
re 
i 


ue 
e, 


ay 


114 


approximation to the. high temperature férromagnetic suscepti- 
bility of the. Ising-model. For example, the asymptotic 
behaviour: of self-avoiding walks can be represented by the 


form 
Cla) (ish) ; (6-2) 


where wu (the "critical point") is called the attrition 
parameter and for.a very long walk represents the number of 
ways the walk can. be continued for one more step. A 
rigorous proof of. the ere. of the attrition parameter 
was given by Hammersley (1957). The numerical values found 
for the attrition parameter wu and the index y in (6-2) 
differ only by a few per cent from the values of the 
corresponding. quantities for the high temperature ferro- 


Magnetic susceptibility of the Ising model. 


| a: 
-~lLiqeckus oitenpsmoxsa2 meersar 


npsorqmyan ond ,alqiaxge Wyott | tobe 2m. 
or yd KesneRerges od ‘hao extew ‘tite 


‘Cm~a) : ~ (my~£) 689, 
: ff, i” ft 4 pres. ve vas 4 ; 
a, ‘ 


to zodaudi ate sidensugey dLaw ena cae eee 

ee able one gete, “GR hawn.tainoe fed Ans: 3 i 

| toJemeteq sobs izsds, eat ee sonetelxe: re 
biratat lav isctuena oem. . (seen) ane 
(S-o) Ad Y xebat eft ae i ‘steno nods $4 

ad? Yeo seuvlev | ons most gaso 360 wo a 

“pxne? Siteexegmet riehd odd, 10%, pidicie 

»Lehou. eadel ont. a 


~ 
a 
. N 
a 
Ww ii 
4 Fi 
i oh f 
: h 
. 4) ’ 
Pal 6 - 
‘s A 
Pe \ 
7 i | ; 
ye 
a 


nad wank ha 


rt 
, 
* ’ 
=. A " 


115 


B. Enumeration of Self-Avoiding Walks 


The enumeration of self-avoiding walks is usually 
carried out in one of two ways, a sampling of walks of very 
long length via Monte-Carlo techniques (Gans 1965) or an 
exact enumeration: of the first few walks (Martin et al. 1967). 
In either case the problem can only conveniently be attacked 


with the aid of a computer. 


We have carried out an exact enumeration of 
self-avoiding walks on the hydrogen peroxide and hyper- 
triangular lattices. Such a calculation is most easily carried 
out by means of a chain counting theorem derived by Sykes 
(1961). In this: section we shall briefly outline the 
derivation of this theorem along the lines originally 


followed by Sykes. 


We focus our attention on a self-avoiding walk of 
(n-1) steps on a lattice of coordination number q. If we 
add one more: step. to. the walk in any one of the (q-1) =,0 
allowable directions, one of three things can.occur. We. 


can form a self-avoiding walk of n._ steps, or one of two 


topologically distinct linear graphs (Fig. 6.1). 
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Bs Sy 


Qt+b=Z 
(a) (b) (c) 


Fig. 6.1. The result of adding a step to a self avoiding 
walk of (n-1) steps. (a) Self-avoiding walk of x2 steps. 


(b) Tadpole, Ue: (c) Polygon, Po: 


In Fig.: 6.1, graph (b), which:is known as a 
tadpole, results from a self intersection a steps from the 
origin. The number of such graphs on a lattice is denoted 


by.g sr" 6 Every. tadpole will be walked twice, the head 


ayb- 
being walked once in each sense. Graph (c) results from 
a self intersection at the origin. The number of such 


graphs per site is denoted by ppg , the number through a 


given point on the lattice is given by LD > . The & steps 
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forming the’ polygon: will also be walked twice, once in each 
sense. We can combine the above statements into an equation 


of the form 


which gives us a simple recursive relation between c and 
ey “in terms of tadpoles and polygons of order & 
We can. extend the idea of the previous. paragraph 
and use the same technique to count the tadpoles. If we 
add a step to the.tail of a tadpole Tyo. y’ the result will 
lA 


be either a new tadpole. Ty y or one of three topologically 


distinct linear graphs shown in Fig. 6.2. 


(a) (b) (c) 


Fig. 6.2.’ The three topologically distinct linear graphs 
that result. fromthe .addition.of a step to the.tail.of a 
tadpole. (a) Dumbell, (a,b,c) s)- (b) Figure eight, (a,b) 
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(c) Theta graph, (a,b,c) > 
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For a given: lattice the number of dumbells, figure eights, 


and theta graphs.per: site are denoted by (a,b,c) (aD) i, 
8 


db’ 
and (a,b,c), respectively. Thus by eliminating the 
tadpoles the number of self-avoiding walks of order & can 
be made: to depend:.on: four types of topologically distinct 


linear graphs, namely polygons, dumbells, figure eights, 


and theta graphs. 


The mathematical steps involved in eliminating 
the tadpoles from the recursive relation (6-3) are straight 
forward but rather lengthy and we shall not reproduce 


them here, but quote only the final result, i.e. 
Copp 20S gta cy» = 2h p -2( M41) py, +8 phy (BrP?) ab 


ee Ee, ) ie ye (aby cya (6-4) 
8 
M+] L+1 


where & > 1 and the summations are over all dumbells, 
figure eights, and theta graphs of order %+1l on the 
lattice. Equation (6-4) is known as Sykes’ chain counting 
theorem. The great advantage of this theorem is that the 
polygons, dumbells etc. are much less numerous than the 
walks themselves, resulting in a large saving of computer 


time, which enables one. to calculate walks of higher order 
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than would otherwise be possible. 


Equation (6-4) can be cast into the alternative 


form 
C(x) = 1+(1-ox)* [qx(l-ox)-2(1-x) } % p, x° 
r 
+ 8 } (a,b,c) Sela: (a,b) a 
Le wy db } ’ P 
u 
| mee g A929 
tI 2) (ag b ic) an 1. (6-5) 


& 


It is perhaps worth mentioning that (6-5) is not simply a 
restricted form of the susceptibility counting theorem of 
the Ising model (c.f. eqn.(3-31) Chapter III). If one 
restricts the susceptibility in the Ising model to a sum 
over polygons, dumbells, figure eights, and theta graphs, 
the counting weights for polygons and theta graphs differ 
from the counting weights for the same graphs in (6-5). 

In both. cases ,(3-31) and (6-5), the counting weights.for 
dumbells and figure eights remain the same. This simply 
points out the well known fact that although the expansions 
for self-avoiding walks and the Ising model susceptibility 
appear to be similar, there are fundamental differences 
between the two problems, namely that the counting weights 
in (6-5) not only depend on the degree of the vertices 


& 


but also on the topology of the graphs as well. 
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Using (6-5) we have enumerated walks of up to 
thirty steps for the hydrogen peroxide lattice and walks of 
up to fourteen steps for the hypertriangular lattice. The 
self-avoiding walk generating functions for these two 


lattices are given by 


(Hydrogen Peroxide) 

Cea alet Sxet 6x- 0+ Lex? 24x" s+ 48x> 4° 96x° +:°192x7 
+ 384x° + 768x?°+ 1506x*® + 2982x'! + 5904x?? 
+ 11688x*? + 23094x’* + 45678x** + 90000x!® 
+ 177660x'’ + 349938x!® + 690192x'* + 1359288x?° 
+ 2678808x*). + 5271558x7* +.10381926x"° 
+ 20419224x7* + 40191204x** + 79025742x** 
+ 155470668x*’ + 305587564x** + 600950160x*? 


+ 1180825386x°?" +... (6-6) 
and 


(Hypertriangular) 
C(x) = 1 + 6x + 30x* + 144x? + 696x* +3330x° + 15774x® 
+ 74484x’ + 351192x® + 1651806x” +. 7753182x°*° 
+ 36342882x!* + 170163366x'* + 795893904x* ’ 


+ 3729252820" * + 234 =. (6-7) 
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The number of polygons, dumbells, figure eights, 
and theta graphs of order & embeddable in the hydrogen 
peroxide and hypertriangular lattices are listed in 


Appendix J. 
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CHAPTER VII 


ANALYSIS OF SERIES EXPANSIONS 


A. Methods of Analysis 


The thermodynamic functions of the three dimen- 
sional Ising model, initial susceptibility, specific heat, 
Spontaneous magnetization, etc., are known to us only 
through their power series expansions, which for the 
hydrogen peroxide and hypertriangular lattices have been 
derived by methods outlined in previous chapters. From 
these series we hope to locate the position of the critical 
point and svedi et the behaviour of the various thermodynamic 


functions near the critical point. 


Power series expansions for the Ising model 
generally fall into three categories (see C. Domb in 
Proceedings of the International Conference on Phenomena 
near Critical Points 1965): 

(1) The coefficients are all of the same sign, 

(2) The coefficients alternate in sign, 

(3) The magnitudes and signs of the coefficients 

exhibit more complicated behavior than (1) or 
(2). 


In (1) the dominant singularity (the physical singularity or 
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critical point) lies on the positive real axis, while in 

(2) the dominant singularity lies on the negative real axis 
and in (3) the dominant singularity lies elsewhere in the 
complex plane. We are primarily interested in series of type 


CL )isior, (2) ¢ 


For the Ising model of a ferromagnet in zero field 
we expect. to find a singularity on the positive real axis 
(the Curie point or critical point) and possibly a singularity 
on the negative real axis corresponding to the antiferromag- 


netic Néel point. 


The power series for the high temperature initial 
susceptibility of the spin one-half Ising model ferromagnet 
is of type (1) for all regular lattices and experience has 
shown this series to be the most useful for obtaining esti- 
mates of the critical point. The power series for the 
spontaneous magnetization of the hydrogen peroxide lattice 
is also of type (1), the only three dimensional lattice other 


than diamond for which this is so. 


Assuming that the asymptotic form of the thermo- 
dynamic functions of the three dimensional Ising model is 
the same as that for the two dimensional Ising model (see 
hap . Le 4d kes 
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there are two methods which have been shown capable of 
giving accurate results from an analysis of the power series 
expansions. These two methods are the ratio method (Domb 


and Sykes 1961) and the Pade approximant method (Baker 1961). 


The ratio method follows directly from (7-1), 


where the ratio of successive terms is given by 


where lg = 1/Xe. Provided the assumption (7-1) is valid and 
the convergence of the series is fairly rapid, a plot of 
an/an-] against L/n should yield a straight line. The inter- 
cept of this line with the (l/n) = 0 axis gives the critical 
point ug, and from the slope one can find the critical index 


Pe Le Co slope = Uc(p-1) © 


Given an estimate of the critical point we can 
ebtain successive estimates of the critical index p from the 


relation 
(p), = Ltn [(up/ug) - 1, (7-3) 


where up, = ap/an-]- If we have an accurate estimate of the 
critical index p, a refined estimate of the critical point 


can be obtained by calculating the successive estimates 


(Ugly = (nun) /(n + pl). (7-4) 
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Provided we have accurate estimates of both the critical 
point Uo and the critical index p we can obtain successive 
estimates of the amplitude A by factoring out of the nth term 
of the power series expansion that part arising from the 


binomial expansion, i. e. 


tie enileane a) de: (7-5) 


The Padé approximant method, which was first 
applied to the Ising model by Baker (1961), approximates a 
function by the ratio of two polynomials. The (L,M) Padé 
approximant to a power series is defined by 


W (x) 4 QL (x) = bo+tb1xtb2x2 ++ 0. ¢by,x" 
Pm (x) CotCpXtC2x® +o ° °+CyxX 


0 (7-6) 


The coefficients b,, bl.,°°°, by, Cor Clr*** sy Sy are calcu- 
lated by requiring the expansion of the right hand side of 
(7-6) to agree exactly with the given power series, W(x), 
through order (L + M) < R, where R is the order of the term 


at which W(x) terminates. 


In order to use the Padé approximant method most 
effectively for the type of functions we are dealing with 


one takes the logarithmic derivative of (7-1), 


(a/dx) ln W(x) = =p/(x-X@) , (7-7) 
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a process which converts the singularity into a simple pole. 
The (L,M) Padé approximant has M simple poles in its denomin- 
ator and hence there is the possibility that the Padé 


approximant may be able to locate the pole we are seeking. 


If we have an accurate estimate of the critical 
index p we can obtain estimates of the critical point by 


forming Padé approximants to 
1 
w(x) ]7/P » -al/P x /(x-xQ) (7-8) 


which has a simple pole at x = Xoo 


Having obtained estimates of the critical point 
from Padé approximants to (7-7) or (7-8) we can use that 
estimate to obtain estimates of the critical index p by 


forming Padé approximants to 
(X-Xo) (d/dy) in W(x) = -p (7-9) 


and avaluating them at the critical point x = X,. Another 
technique for obtaining estimates of the critical index was 
introduced by Baker et al. (1967). This method consists of 


forming Padé approximants to the function 
(d/dx) in [(d/dy) W(x)]/(d/dx) In W(x) = (p + 1)/p (7-10) 


and evaluating them at the critical point. 
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One can also use the Padé approximant method to 
calculate estimates of the amplitude. Taking Padé approxi- 


mants to the function 


(x=%—) (W(x) ]27P we -al/P x (7-11) 


fon 


and evaluating them at the critical point we can get 


estimates of A. 
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Bed The Ising Model 


The series expansions derived in Chapter III and 
IV have been analyzed using the methods outlined in the 
previous Beation’ The results of a Padé approximant 
analysis are most conveniently displayed in the form of a 
table consisting of a rectangular array of the (L,M) 
approximants. As we have several such tables they have 
been grouped together to form Appendix K. The results 
of the ratio analysis are displayed in Fig. 7.1 and 


Table 7.1. 


The series expansion for the high temperature 
initial susceptibility of the hypertriangular lattice was 
analyzed first. Since it has Been rather well established © 
for some time that the susceptibility critical index 
y = 5/4, we have obtained estimates of the critical point, 
Vo, from Padé approximants to ened we as shown in 
Table K.1. Without assuming a value for y, estimates of 
Ve were pine Hae from (d/dv) 1ln x(v), and are listed in 
Table K.2. The two methods yield results consistent to 
five decimal places. Taking as our best estimate of the 
critical point of the hypertriangular lattice, 


Vo = 0.222087 + 0.000005 from Table K.1, estimates of 7 


were obtained from Padé approximants to (ve-v) (d/dv)in x(v) 
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and are listed in Table K.3. Estimates of y obtained from 
Padé approximants to (d/dv) [ln (d/dv) x(v)J/(a@/av) In x (v) 
are listed in Table K.4,. The results of the two methods 
are guite consistent and to four decimal places the entries 


of Table K.3 are consistent with y = 5/4. 


The amplitude, Cy (v<), of the susceptibility 
Singularity in x(v) = Cy (1-v/vg) * has been estimated 
from Pade approximants to (Vc-¥) [x(v) ] 475; the results 
are shown in Table K.5, Pade approximant estimates of the 
dimensionless free energy per site and internal energy, 
evaluated at the critical point, are given in Tables K.6 
and K.7 respectively. Finally, estimates of the amplitude, 
Ay (Ke), of the specific heat singularity in C (K) # A, 
(1-K/K,) — were obtained from Pade approximants to 


L/ (2+a) 


(K-Ke) {(a?/dK?) C (K)} and are listed in Table K.8. 


The high temperature series expansions for the 
hydrogen peroxide lattice have been analyzed in the same 
way as those for the hypertriangular lattice, with one 
important difference. Our best estimate of the critical 
point ve of the hydrogen peroxide lattice, taken from 
Table K.9 is ve = 0.51815 + 0.00006. However, using our 
best estimate of the critical point of the hypertriangular 


lattice, ve = 0.222087 + 0.000005, and the star-triangle 
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relation (5-12), we get v, = 0.518140 + 0.000006 for 

the critical point of the hydrogen peroxide lattice. It is 
this latter estimate that we have used in calculations 
requiring an estimate of the critical point of the hydrogen 
peroxide lattice. As an indication of the precision that 
has been achieved in estimating the critical points of the 
hypertriangular and hydrogen peroxide lattices, it may be 
noted that the best estimates obtained from Tables K.1 and 
K.9 satisfy the relation (5-12) to five figures with an 


error of only 4 parts in 10°. 


We have also used the ratio method to analyze 
the susceptibility and specific heat series of the hyper- 
triangular lattice and the susceptibility series of the 
hydrogen peroxide lattice. A plot of bn/bn-] against 1l/n 
for the susceptibility of the hypertriangular lattice is 
given in Fig. 7.1, where g = y-1 and pee l/vco-. A com 
parison of the estimates listed in Table 7.1 with the 
corresponding estimates obtained from Pade analysis shows 


that the two methods have yielded consistent results. 


For the hydrogen peroxide lattice it was not 
possible to obtain any meaningful estimates for the 
specific heat singularity amplitude. In principle it is 


possible to obtain ve from aad using (5-46). However, 
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Fig. 7.1. Plot of successive ratios b,/b,_} 
against 1/n for the initial susceptibility 
of the hypertriangular lattice. te= 1/ve, 


g = y-l. 
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we have only rough estimates of aut (c.£. Tables K.8 and 


7.14) and so we can do no better for ae 


Precise estimates of the amplitude of the 
susceptibility have been obtained for both lattices and 
thus a comparison with the ratio predicted by (5-51) is 
possible. From (5-51), taking te = 0.225850 and 


KEP = 0.573794 we find 


HP HT = 
Cl" (Ke) /Cy* (Ke) = 1.248055, 


whereas from separate estimates from Tables K.5 and 


K.13 we find 
COP OAK JI/Cr AR) ae Ae 2488s 


where 


Cy (Ko) = (sinh 2Ke/2Ke)’ Cy (vg). 


The coefficients in the low temperature series 
expansions for the hydrogen peroxide lattice are all of 
the same sign and thus the dominant singularity is the 
physical singularity or critical point. As pointed out in 
Section A of this chapter the ratio method is applicable 
to such series, However, the convergence of the various 
low temperature series expansions for the hydrogen peroxide 


lattice has proven to be painfully slow and as a consequence 
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it has been impossible to draw any conclusion about the 
behaviour of these functions from a ratio analysis. By 
necessity then we have been forced to rely solely on the 


Pade approximant method. 


In Table K.16 we have listed for the hydrogen 
peroxide lattice, estimates of the critical point 
Zo = exp (-23/kTg) obtained from Pade approximants to 
(d/dz) ln I(z). A critical point given by Ze = 0.316 
+ 002 is consistent with the entries of Table K,16 
However, our best estimate of the critical point of the 
hydrogen peroxide lattice is given by ve = 0.518140 
+ 0.000006, or since ze = (l-ve)/(ltve) by zc = 0.317402 
+ 0.000005. It is this value that we have used in all 


calculations requiring an estimate of Zc. 


Estimates of the critical indices a', Band y', 
the critical indices of the specific heat, spontaneous 
magnetization, and susceptibility respectively, have 
been obtained by employing the Pade approximant method 
as given by (7-10). The estimates of 8 given in Table K.1? 
are consistent with the value § = 0.305 + 0.005. The 
estimates of y' listed in Table K.18 show a fair amount 
of scatter and no one particular value occurs with any 


frequency, However, as it is generally thought that the 
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diagonal and near diagonal Pade approximants are the most 
reliable, we estimate y' = 1.27 + 0.02. Estimates of a! 
given in Table .13 like those of 8 are fairly consistent, 


and we take as our best estimate a' = 0,23 + 0.01. 


In Table K,20 we have listed estimates of 
(~£p/KT.) (the dimensionless free energy per site) obtained 
from Pade approximants to £/kT evaluated at the critical 
point. As our best estimate we take (f,/kT,) = 0.931, 
which to three figures agrees exactly with the estimate 


obtained from the high temperature series. 


We have also calculated estimates of the critical 
energy by forming Pade approximants to U(z)/J and 
evaluating them at the critical point. These results 
are shown in Table K.21, The entries in Table K. 21 seem 
to indicate a value for U,/kTg of = 0.51. However, it 
should be noted that the (8,8) Pade approximant,a diagonal 
approximant making use of all the terms of the series, 
is in good agreement with the high temperature result 
U-/kT, = -0,.479. Since the entropy is related to the free 
energy and the internal energy through (3-28), our best 
estimate of the critical entropy for the hydrogen peroxide 
lattice obtained from low temperature series expansions is 


S./k = 0.458, as opposed to the high temperature estimate 


of So/k = 0.451. 
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Cc. The Self-Avoiding Walk Problem 


The generating functions for self avoiding walks on 
the hypertriangular and hydrogen peroxide lattices have also 
been analyzed by the ratio and Padé approximant methods. In 
Fig. 7.2 we have plotted the ratios of successive terms of 
C(x) for the hypertriangular lattice. Nearly all the points 
lie on a straigth line, which indicates that (7-1) is pro- 
bably a valid assumption for this problem also. From 
Fig. 7.2 we estimate that »p = 4.6179 and g = 0.1665, from 
which we conjecture that g = 1/6 exactly, in agreement with 
estimates from other three dimensional lattices (Essam and 


Sykes 1963, Martin et al. 1967, Guttman et al. 1968). 


In Table 7.2 we have listed estimates of the 
attrition parameter yw, ana the singularity amplitude A, for 
the hypertriangular and hydrogen peroxide lattices. The 
estimates of the attrition parameter are calculated from 
(7-4) and the estimates of the amplitude from (7-5). In 


calculating both u and A, we have assumed g = 1/6 


n 
exactly. For the hypertriangular lattice the limiting value 
of uw in Table 7.2 agrees very well with the value of u 

taken from Fig. 7.2, which indicates that the conjecture 


g = 1/6 is probably correct. 
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Fig. 7.2. Plot of successive ratios Cy/cy.} 
against 1l/n for the self-avoiding walk 
generating function of the hypertriangular 


lattice. 
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TABLE 7.2 


Ratio estimates of the attrition parameter 


Hypertriangular 


Hydrogen Peroxide 


n Un An n 5 Cuntin-a) An 
7 4.612135 1.09898 20 1.95417 1.1226 
8 4.618775 1.09894 a7 1.95419 1.1216 
9 4.617910 1.09870 22 1.95415 1.1194 
i LO 4.616814 1.09819 23 1.95416 1.1184 
11 4.617517 1.09786 24 1.95424 1.1163 
12 4.618026 1.09764 25 1.95426 1.1153 
13 4.618029 1.09742 26 1.95450 1.1135 
14 4.618074 1.09722 27 1.95450 121125 
' ~~2 --- 28 1.95460 1.1108 
‘ a --- 29 1.95462 1.1099 
' ort “<> 30 1.95469 1.1083 

oe 4.6181 + 1.095 + ' “a= “ao 

0.0046 pyarey ta 

a on om an ale eos 

§ eae ane taxe ee exo come 
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0.020 Oat 
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ratios of successive terms of C(x) for the hydrogen peroxide 
lattice, we have calculated an average of successive pairs of 
estimates of u. For the same reason we have not attempted to 
estimate yu and g for the hydrogen peroxide lattice from a 


ratio plot such as Fig. 7.2. 


We have also obtained estimates of the inverse of 
the attrition parameter using the Padé approximant method 
(see TableK.22 and Table K.23 of Appendix K). These estimates 
are quite consistent with those obtained from the ratio 


method, 


On the basis of successive estimates of the attri- 
tion parameter and singularity amplitude shown in Table 7.2 
our best estimates for the hydrogen peroxide lattice are 
u = 1.956 + 0.020, A= 1.10 + 0.11, and for the hypertri- 
angular lattice they are yp = 4.6181 + 0.0046 and A = 1.095 
+ 0.011. 


f 7 yr a. bait ie, ee zm 


5 


i ; “it ae 
ob iuesbq neyporbyd @tis aa8 wa s 
) wees ai: er 


tn sitaq avireroouet Yo omnsovs ape! 


” 


ot) botquadts tod evad ov none 


by f 
"9 aenornd ef 2o, setanigem amtegdo moa 
husidom taamizosgap Shee. wy pentew 
eotenicas geerit | «idea * thn oldu x ban § 
otoes eld dost bentetde sand dike steianoe @4 
wake it & ne ne a 


haw 
~ bases ott Yo aataanides sninpioatae ie etesd ¢ oils oO as 
an Ab Gy 


4 vit pr 4 iu, tM) ye , 

Cif afdet at awots aha hha gure mR mies a es 
fi t, Set 4 ah : an WT ” 

oun edhe sat sbkccxeg moposbyd ef 0% se : ao 38 


wbad ened oft 20% Gne hte se ‘Ohet « e A 


ole hans : 


20.4 » A bate atin. A ee ee 


al wd A 


a af Wind y 7 py) 


Brey mney Pet Ys 
aie La ah 7 i a 
v4 ' eave. F ye 7 “ rv ka | 


i : a oe : Ki; iA ui ™ : rs " 
b H ar } : + : y, hi ., } i} 7 
} . ‘ rin as aa , M Y ; 
I F : aa ts ; : 
| ‘ on ! ey i Nt 
ae No A Ae oc | ; 
Ai ) : j ; ; ¥ | - ps bora 
1, of Ai nn ary 
; aes ii Woes Ta Pah Paty, a j f : 
i‘ 1 fi ak P ne » 7 sm 


140 


CHAPTER VIII 


REVIEW AND DISCUSSION 


A. The Ising Model 


Our primary goal in this work as stated in the 
iscrocaction was to obtain improved estimates of the low 
temperature critical indices of the three dimensional Ising 
model. It is the opinion of the author that we have obtained 
improved estimates of the low temperature critical indices, 
though admittedly we have not obtained the precision we 
had hoped for. However, we intend to pursue this problem 
in the future with the hope of extending the low temperature 
series. This combined with a more sophisticated analysis 


should shed further light on this problem. 


The first stage of this work consisted in deriving 
high temperature series expansions for the hydrogen peroxide 
and hypertriangular lattices. The reason for this was to 
obtain a highly precise value for the critical temperature 
of the hydrogen peroxide lattice. This precision was 
achieved by first obtaining a highly precise estimate of the 
critical temperature of the hypertriangular lattice from the 
series expansion for the initial susceptibility and then using 


the star-triangle transformation to obtain the final estimate 


a. | ae ouce, te i - rev ; un POM bid ” 4 Ry iid ans ae ad 
’ oe a : Ue ha vee ™ ; 


oor ) vr | . hi bi 1) 
ran MRA, | 


oe ¥ (* Pack i wor on 


7 ne - ey) Tar 


~ ~ 


ot nt Bedats ee snow atts at Isop. sala 
) 


VREEM, Sai-cttyy Wie 


WO eit Zo wick anilihew bewonget ero ant eew a 

eit it shit cia SiN LP oon.itd ant tom epital Lp is ae 's 

, hal ‘righ 

Hontnads eved ow datz sorters odd ‘0 motaigo erty at 4 

,geoinat Ipsos tvs eres seqmet ‘onmed va | fi 

aw nobe Loene erly begtatde tor evati ow yl + a 
epalcler tes wha oc Rag ad, baezas ow seat 03 

erp weol echt) Qatierr%e 20 eer “es 

a ie yl te botsotsatiqas rat & awe 


iting it i hee 


a 
i) Ey: f 


ny am 


it eM ER Sea wish Wee 

idle kalal nk Bededtanee Yaow eins to opsta tex ri aa 
sbixezbd neyorbyd edd 293 snohansaxs eine as 
ot wow eiad nol soars od? sBeoldseL % an er api 
owiereqeet Leoliiss edd x03 aulev atoexg Hers 
he: my 

ey nosaiodng aldv aokiaat sbixoxeg nop ‘ Bn ye i 

ott to etamiteo watoang winks rt clad 2a ue 


141 


for the critical temperature of the hydrogen peroxide lattice. 


A further result of this work is that critical 
parameters for the three dimensional Ising model are now 
known for the very lowest coordination number, i.e. q=3, 
hence we are in a position to examine the dependence of these 
parameters on coordination number. In Table 8.1 we have 
tabulated most of the known results for the two dimensional 
Ising model. For the three dimensional Ising model we have 
listed the critical parameters in Table 8.2. For a given 
dimension the linear chain represents the lower limit of 
coordination number (q=2) while the mean field theory repre- 
sents the upper limit (g=~). Although Tables 8.1 and 8.2 
represent almost all that is known about high temperature 
lattice dependent critical parameters, the dependence of 
these parameters on coordination number is better appreciated 


when displayed graphically. 


In Fig. 8.1 we have displayed the dependence of 
Vc = tanh(J/kTc) on coordination number q for two and three 
dimensional lattices. The letters labelling the calculated 
critical temperatures are obvious abbreviations for the 
lattices listed in Tables 8.1 and 8.2. The curves are not 
a best fit to all the calculated points but are obtained by 


a simple procedure. They represent the function 
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Fig. 8.1. The dependence of the critical 
nee Vo = tanh(J/kT,), on coordination 
number q for two and three dimensional Ising 
models. Circles indicate calculated points 
for loose-packed lattices; triangles, close- 
packed lattices. The dashed curve repre- 


ees i 


sents the function v, = (q-1) : 
-0.950 


solid curve the function vg = (q-1) 


144 


> 2 


i 
‘ 
, rT 
i 
i 
j i 
-_ it 
[ 6 
, 
7A, 
DD hd 


Kako es 


~seodp ssebgniaind & 


‘ iY 


an y ea 


awe eHs Y 


A a ae 


i i | 
ae +4 any yiats 


a 


it ' 


~ 0.806 
a My (S20) 


~ 0.950 
Ve = (4-17 98 


BCC (1,2) 
$C(1,2),FCc(4,2) 


1.0 


tanh(Ivkk) 


vox 


e.0 
8.0} 
0 


@.0} 


ST aaa 
aL y Ari's dein ae 
Reh TS a Burgi to ee 


. ae we Y eae anh | | 
; #0 AK, a 4 - at ir 
ds ea Raat Aad 


2. 


ea soup(ane £) 


. 7 7 aa 
na Yf cota oa hs Asam ay " 4 
&0 , vi Pe OG ate aa nae { A iy 
i] 


a | : } ~ } 

an ty 

- P an 0, he 4 i = 
+e re : ‘, b “phim sits we ra ‘gyi ‘esenie 7 i Z| 


eos!) pen 2) rey 


me 


| ) | . bts s 
; ’ : ow ne oe | 

if i 5 ; e ff ; 7" 3 a J my F a iN 
"eee Y ne 2) eee | io Lan "an 


cea age — eae 7 a8 eco “ ‘ 


145 


vo(a) = (q-1) “4, (8-1) 


where ag is chosen so that ve(2d) = wae the critical point 


of the simple hypercubic lattice. 


In two dimensions the simple hypercubic lattice 
is the square lattice and a, = 0.806; in three dimensions 
the simple hypercubic lattice is the simple cubic lattice and 
a3 = 0.950. We see that the curve for d=3 is an excellent 
fit to most of the critical points over the whole range of 
q and thus (8-1) has predictive value for higher dimensions. 
Fisher pet Gaunt (1964) have computed ve for simple hyper- 
cubic lattices only in four, five, and six dimensions. 
From their data we obtain a, = 0.984, a, = 0.996, and 


a. = 1.000. Since the completion of this work, Moore (1969) 


6 
has obtained the critical points of the four dimensional 
hyper-simple cubic (hsc), hyper-body-centered cubic (hbcc), 
and hyper-face-centered cubic (hfcc) Ising lattices. Moore 
found viS¢ = 1/6.725 = 0.149, vg°°° = 1/14.510 = 0.0690, 


0.0455. The coordination numbers of 


i 


and vafcc = 1/21.984 
the hsc and hfcc lattices in four dimensions are eight and 
twenty~four respectively. For the hbcc lattice Moore has 

considered only the sixteen nearest neighbours lying on the 


crystal axes, overlooking the eight nearest neighbours lying 


P h hnpcc: — 
off the crystal axes. Hence, taking q chee Pel: as = 16, 
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qhfcc = 24, and using (8-1) we find v8S¢ = 9.149, 


bee 
ve 


agreement with Moore's results and further substantiating 


0.0700, and vitSC = 9.0460, all in excellent 


the usefulness of (8-1) for estimating critical points of 


Ising models. 


The amplitude of the specific heat singularity, 


is rather insensitive to q and also is not well known 

for most lattices. C,, on the other hand, is well deter- 
mined for most lattices and varies appreciably from lattice 
to lattice. However, the variation depends strongly on other 
aspects of lattice weducHice besides q. Hence we have not 
displayed graphically the q dependence of either A, or C,. 
The critical energy and entropy both vary rather smoothly 
with q as illustrated in Figs. 8.2 and 8.3, but no such 


simple expression as (8-1) has been found for these quantites. 


As is well known, the Ising model serves as a 
model for a lattice gas (Yang and Lee 1952) so we have com- 


puted the critical ratios (Fisher 1964(b)). 
Poo/P¢kTo = 2[S¢/k - Up/kT. - q/2K,] - era? 


The dependence of the critical ratio on coordination number 


is displayed in Fig. 8.4. 
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Fig. 8.2. The dependence of the critical energy, Ug, on 
coordination number, q, for two and three dimensional 
Ising models. Circles indicate calculated points for 
loose packed lattices; triangles, close packed-lattices. 
The top curve is a best fit for two dimensions; the 


bottom curve for three dimensions. 
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Fig. 8.3. The dependence of the critical entropy, So, 
on coordination number, q, for two and three dimensional 
Ising models. Circles indicate calculated points for 
loose-packed lattices; triangles, close-packed lattices. 
The top curve is best fit for three dimensions; the 


bottom curve for two dimensions. 
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Fig. 8.4. The dependence of the critical ratio, P,/pckT,, 
on coordination number, q, for two and three dimensional 
Ising models. Circles indicate calculated points for 
loose packed lattices; triangles, close-packed lattices. 
The top curve is a best fit for three dimensions; the 


bottom curve for two dimensions. 
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At this point it might be worth noting that for 
argon S./k = 0.43 (Betts and Ditzian 1968), in dase aeae agree= . 
ment with the:value for the hydrogen peroxide lattice. Also, 
Mikolaj and Pings (1968) have recently argued that the 
coordination number of a liquid is meaningful and that q = 3 
for argon at the critical point. However, the value for the 
critical energy for argon corresponds to q = 5. Nevertheless 
we believe this is sufficient Justis ieatien ae the further | 
study of the applicability of the Ising art on the pretaed 


peroxide lattice to simple liquids. 


Briefly the main results of this work are as 


follows: 


1. Three "new" regular lattices have been intro- 
duced for the study of lattice statistical problems. 

2. A large number of weak lattice constants have 
been obtained for the hydrogen peroxide and hypertriangular 
lattices. 

3. A large number of strong lattice constants 
have been obtained for the hypertriangular lattice. 

4, Exact relations have been derived, valid in 
the critical region, relating amplitudes of (a) specific 
heats and (b) susceptibilities for a triple of generalized 


4 o 
honeycomb, triangular, and kagome lattices. 
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2. The critical temperatures for the Ising model 
on all three lattices has been obtained to high precision. 

6. The previously accepted value for y has been 
confirmed. 

7. “Improved" estimates of the low temperature 
critical indices a°, B , and yy’ have eae obtained, which 
indicate the scaling law hypothesis is not valid for the 
three dimensional Ising model. The estimates found for y 
and § are in excellent agreement with the experimental 
results for beta brass, indicating that the Ising model — 
provides a valid description of the order-disorder transi- 
tion in binary alloys. 

8, The susceptibility amplitude, critical energy, 
and critical entropy have been determined to reasonable 
precision on all eneed lattices. 

9. A simple dependence of critical temperature on 
dimensionality and coordination number has been discovered. 

10. The possible relevance of the hydrogen peroxide 
lattice for simple fluids near the critical point has been 


indicated, 
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B. The Self-Avoiding Walk Problem 


We have derived exact series expansions of the 
self-avoiding walk generating functions for the hydrogen 
peroxide and hypertriangular lattices. From a numerical 
analysis of the series coefficients we have obtained esti- 
mates of the attrition parameter, critical index, and 
amplitude for both lattices. The value found for the critical 
index agrees with the value previously found for other three 


dimensional lattices. 


A considerable amount of effort was expended 
searching for the equivalent of the star-triangle trans- 
formation for the self-avoiding walk problem, since such a 
Beets formation would provide us with a rigorous check of the 
numerical analysis and would presumably also show us how to 
transform the generating function from the hydrogen peroxide 
lattice to the hypertriangular lattice. However, these 
efforts have so far been fruitless. It has been pointed out 
to us by Sykes (Sykes, private communication) that the 
attrition parameters of the hydrogen peroxide and hypertri- 


angular lattices should obey the following inequality, 


v>y u?/(l + p) (8-3) 
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where v and up denote the attrition cevatieteve of the hydrogen 
peroxide and hypertriangular lattices respectively. Taking 

wp = 4.6181 we find v > 1.948 which it just is at 1.956, a 
difference of only 0.4%. The equality in (8-3) holds for 

the attrition parameters of Bethe lattices of coordination 


numbers six and three respectively. 


Part of the work reported on in this thesis has 
also appeared in the physics literature. The high temper- 
ature behaviour of the Ising model on the hydrogen peroxide, 
hypertriangular, and hyperkagome lattices is discussed in 
Leu, Betts, and Elliott (1969). A discussion of the self- 
avoiding walk problem on the hydrogen peroxide and hyper- 


triangular lattices appears in Leu (1969). 
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APPENDIX A 


The first twelve coefficients im the series expansion 


for the zero field partition function of the hypertriangular 


lattice, 

P(1,N) = 0 

P(2,N) = 0 

P(3,N) = N(A) 
= N 

P(4,N) = 0 


P(5,N) =N (QQ) 


= 3N 
mere N)) = Nh Cy te NDS) pe (A A) 
= 1/2 N2 + 14 1/2 N_ | 
P(7,N) = N «DS 
= 33N 
P(8,N) =N ag ) +N (OX ) +N (O Ls:) 
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APPENDIX B 


Weak lattice constants of non-magnetic graphs required to 


calculate the partition function and initial susceptibility 


of the hypertriangular lattice. 


Graph 


Lattice 


Constant. 
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33 
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3 N=-33 


1/6 N2-3 1/2 N + 20 1/3 
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| Lattice 

Graph Constant 
(xX ) 51 
(CYA) 33 N-462 
(Od) 4 1/2 N-76 1/2 
cas 3846 
(00 ) 459 
(CH ) 651 
(QA4 ) 60 
(Q) 15 
(POM ) 15 
(CO ) 18 
(ja) 72 N-1134 
(Od) 45 N-840 
(() Mt) 9 N-135 
(HA) 15 N-210 


(O4A)11/2M - 19 1/2N + 300 
7e LN 


ae 13229 
(> ) 2664 


(Cx ) 972 


ve 
wae 


enlisted 
FA fKs BOSD 
a4 


Lda LE 


exe a i a emu Me ie \ 


&! 

Bs 

bELiinenm SX 
ChG~u ab 
Mi-a & 
gikew 24 


ehser 
‘boas 
ere 


Lattice 


Constant 


1179 
408 
120 


293 N-5079 
105 N=1668 
N-15 


12 N-180 


157 


Lattice 

Graph Constant 

(OO) 99 N-2037 

(OO) 112 1/2°N-2353 1/2 

(Cp) 45 N-804 
(DD) 41/2 N-76 1/2 


(QAA) 7-1/2 N2-232 1/2 N 
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Weak lattice constants of magnetic ‘graphs which are required 


to calculate the initial susceptibility of the hypertriangular 


lattice from Sykes' theorem. 


Graph 
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APPENDIX C 


Symbolic equations for calculating the lattice constants 


of separated non-magnetic graphs for the hypertriangular 


lattice. 

N(A) N(A) = N(A) + 2N (AA) + an (AA) 
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Symbolic equations for calculating the lattice constants 


of separated magnetic graphs for the hypertriangular 


lattice. 
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APPENDIX D 


Weak lattice constants of polygon, theta, alpha, beta, 
and gamma star graphs for the hydrogen peroxide and hyper- 


triangular lattices. 
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Weak lattice constants of p graphs or polygons 


ann nr eS ASS SRS SC Ss SSS Sn SSS SSA Ss SS SS SSS SO ee sn aerenoemety veseeeeerdascecemeeweasees 
A SR SSS ss SNe rn ese een —eennss 


hydrogen peroxide lattice hypertrianguiar lattice 
g | (2), g (2) 5 
10" 11/2 3 1 
je |S 0 4 ) 
14 21/2 5 3 
16 10 1/2 6 15 
18 23 1/2 7 33 
20 45 8 72 
22 1781/2 9 293 
24 5871/2 10 1158 
26 1579 1/2 11 3846 
28 5026 1/2 12 13229 
30 16,603 1/2 13 50145 
14 192054 


*The first p graph embeddable on the hydrogen peroxide 
lattice is the decagon and no graphs of odd order are 


embeddable on it. 
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Weak lattice constants of @ graphs 


for the hydrogen peroxide lattice” 
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a b c L (a,b,c) g a b c £ (a,b,c)¢ 


2 8 18 28 1047 2 8 20 30 4458 

2 12 14 28 252 oie) ero 30) 522 

4 6 18 28 990 2 14 #14 30 948 

4 10 14 28 150 4 6 20 30 3300 

e Le, k2- 28 237 & 0 16 30 234 

6 8 14 28 96 4 12 14 30 870 

6 10 12 28 222 6 8 16 30 390 

8 8 12 28 78 6 10 14 30 564 

8 ZO x0 , 28 78 6 12 12 30 342 
8 8 14 30 249 

1 9 09 «629 2028 8 10 12 30 366 

lL 23 i> 29 489 10 10- 10 30 42 

3 7 Lg: 29 1890 

3 it iS 2g 230 

3 g3 i3 29 786 

5 5 19 29 984 

5 9 25 29 144 

a dl 13 29 531 

7 7 15: 29 48 

7 9 13 29 216 

7 yi dL 29 135 

9 9 se 72 


“Many 8 graphs which are by inspection obviously not enibed- 


dable in the hydrogen peroxide lattice are not listed. 
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Weak lattice constants of 6 graphs for the 


hypertriangualr lattice 


m2 % 7 15 1/i3@ oie 7392 
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*Many 6 graphs which are by inspection obviously not embed- 


dable in the hypertriangular lattice are not listed. 
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Weak lattice constants of a graphs for the hydrogen 


peroxide and hypertriangular lattices 


a b “peseoeacaned: Umea SARUM “MRAM SUNRN F422 tec T ACI EY 


(i) hydrogen peroxide lattice 
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L 2 l 6 7 8 25 24 
L 6 1 6 3 8 25 12 
L 6 2 7 p. 7 25 24 
2 3 3 8 4 5 25 6 
1 1 1 7 8 8 26 36 
L q Ll 7 2 8 26 24 
y 2 2 8 6 6 26 3 
1 8 Ll 8 pe 8 27 13 
(ii) hypertriangular lattice 
L 2 L 2 2 1 9 2 
Ll 2 L 2 3 1 LO 12 
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Weak lattice constants of 8 graphs for the hydrogen 


‘peroxide and hypertriangular lattices 


a b C a Bites £ g (a,b; 


c,d; e,£) B 
(i) hydrogen peroxide lattice 
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(ii) hypertriangular lattice 


a b Cc d e £ £ (a,b; c,d; e,f) 2 


3 2 L 2 1 2 9 L5 
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L 4 1 4 i 2 aL V2 
d. 2 2 3 l 2 ll 36 
1 2 2 4 1 d ll 36 
| 2 3 3 1 1 a 24 
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Weak lattice constants of y graphs for the 


hypertriangular lattice 
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The first nine coefficients, calculated by the direct method, 


in the power series expansion of the high temperature 


ferromagnetic initial susceptibility of the hypertriangular 


lattice. 
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Lattice constants per site of weak graphs required in the 


calculation, by the direct. method, of the high temperature 


ferromagnetic initial susceptibility of the hypertriangular 


lattice. 
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Lattice Constant 
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Sympolic equations used to calculate the weak lattice 


constants of separated configurations required in the 


calculation by the direct method of the high temperature 


ferromagnetic initial susceptibility of the hypertriangular 


lattice. 
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APPENDIX F 


first seven low temperature ferromagnetic polynomials 


the hypertriangular lattice 
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Low temperature ferromagnetic series expansions for the 


Ising model on the hypertriangular lattice. 
free energy 
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APPENDIX G 


Partial Generating Functions for the Hydrogen Peroxide Lattice. 
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APPENDIX I 


Low Temperature Ferromagnetic Polynomials L_‘2Z) for the 


Hydrogen Peroxide Lattice. 
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APPENDIX J 


The number of polygons (p graphs), theta graphs, dumbeii 


graphs, and figure eight graphs on the hypertriangular 


lattice, 
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Weak lattice constants of dumbell graphs for the hyper- 


triangular lattice. 
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Weak lattice constants of figure eight graphs. for the hyper- 


triangular lattice. 
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APPENDIX K 


Table of Pade approximant estimates of critical parameters 


for the hypertriangular and hydrogen peroxide lattices. 


. P70 i ai ‘ 
ba Ri en | hee 
a, ery. phy ‘omar 4 
' rot 


ao 


sae ae 


244 


GE8OTCC" D 


8680277 °0 


SoTa3eT reTRhuetzazadAy ayy fo 


780CeC-~0 
vVtoccc *0 
I680CCC 0 


BOSCCCC GU 


gh UA)X] 


€L8077Z°0 


ELSUCCC 0 


Cocucce @ 


Stevucec @ 


SScUEcc °C 


B8980CCC -C 


LEsvccc © 


GS30¢CcC- Oo 


SScoccc G6 


LO8027Z°*0 


03 sjZUueUTxordde:spegq wozry pseuTteqgqo 


x 2 


T°M TIAL 


ceeucec 0 


€C680C7C 0 


SSS0ccc= 


"OOSCCCe 


9O060CCT* @ 


EGSOSCE 


0 


“A ~f3uTOd-TeoTtyATrTD syy jo soeyeutasy . 


2 iS oa 
i. 


= “BeaGess.0 ©«=—_- aeBGSSS.0 


S¥SOSS¢ .6 


/ 
= beLosgs.0 


Beoccs. 0 


eesoess .0 
ef 2ocsc = 6 
SOBOSss.0 


TTsosss..0 


vesosss. $ 


{ 


@e80sSs .0 


{rsesss .9 


\ 


PasOsesie 


feeosss +o 


S@ROSSS 0 - 


_ BSeRSS5 .0 


245 


CLVOCCC 0 


SOTAIeT AeTNHuetazzzedAy sya jo 


Seseccc' 0 


WeESOCEC O 


€T60C7C7*0 


OLL0ET? 0 


9ESOTET* O 


Seseccc 0 


OLS0C Ce ~0 


CcOILOcc 0 


V8L0777° 0 


68077272 “0 


GSE8CTCT°O 


9EITECT  O 


ESEECCE~O 


Secelccce o 


ECV8OTCT 0 


CEBSULCC -O 


SVIECCE -0 


SCLECCS 0 


‘VESOCCC 0 


ESSECcc. 0 


EeLicce 


(A) X\UT = OF sjZUueuUTXxOITdde epeq worF pouTeqgqo 


S 
: A 


¢°M = GIEGVL 


SlhICCany 


‘qZutod [TeoTzTrAO |yA JO sejeutjasy” 


bst m 
> £3645; a go greens 


“seaicss-6 
ovegsss,.9  eticss.6 


ceRoets. 0 aSaLlSES 9 -eSvLSss.0 


1 


. 


feesaesk.o scadsse.e aeessse.O SeaNsce.o VEN ESSS.c 


- 


el 


besotce.o 5 peeeses ee errsese.oy | 


agvace?.0 BESO afeesss.d 


246 


T6602 °T 


T667C'T 


SLU aceeb 


v6Eerc’T 


(A)X UT 


AP 


96672 T 


Soore 


C6602 T 


T667c 1 


6L6Vc 0 


L80C77 "0 


P 


€667c°T 


S6677°T 


966772 °T 


VEGVC TL 


CBCSaC it 


A = A ‘32 pezenTerd 


V660C °T 


L6O6vC°T 


v6éEerce Tt 


SI6UV-CET 


TL Ora Caer 


V660c°T 
T667d 1 VLO6VC T 


VVC S.Gae 


(a-"A) 03 squUeUTtTxozdde epeq wory peutezgo 


SoTAReT azetnbuetaqgzeddy oyq go ’ A ‘xepuTt TeoTaTIO oui so soqeUTasY 


e°ua WIdvs 


a ae uP 
»* : 


ay- — 


esses. i 


i - « = = 
ee a = 
be ~ 


beies.t mee A: 
qeebe f geese. f sees .f EROS LL ae * 
u : ; ; ; = nap a 


ae i eeess.i TeORS of a ae 


aa } bee RS . { . = : 


eR, 


7 = <a 


a 7 LL LLL LA LAL ELLA LLL LLL ALLE ELA 


ee om r 7 = 
i : z 
"> eS 


247 


ss SRN SSS SSS SS SS ———————— 


Lvose’T 8 
TE6re Tg CLove’T. Lv0SZ°T L 


9ECGU aE CS6PC'T LE OY. TE 988PC°T O28he°T 9 


9€67C°T 8S67Z°T 9E0S7°T cecse T $96S2°T S 

9S0S2°T €96F2°T 9E607°T O9ES7°T vVCE6Z°T v 
9S0S2°T OC8PC°T BLE9S*T e 

6SS972°T @ 


L80@22Z°0 = PA = A 


; =P AP Ap 
Re pe ,enTeas [ (A)X uT] af tLe —pl Ut} 0} squeutxordde sped wozrzF 


peuteqqo seTa3eT zernbuetazxzedAy eyy go ’ A‘xeput TeoTATIO eyy Jo sezeUTISg, 


fa ATAVG 


ia el => oo yor a ~_ 


nies = ¢ Lopasi egy 
= 


:  MeOSSS. 0 = - : 2 


eazss.f 


OS035.£ ereas,. s 
SERRE. i «=—«-«ETE,L SHEE 
Se0es.i «© - SORES. 


ceeec i 


248 


GLLESO> ET 


9TL6E0 T 


9TLEEO°T -6TLEEOQ T “ELS6tO T 


OTLG6EO°T -GTLG6EO*°T Ze96EO°T —SELE6EO T- Set6c0 = 1 


STICCO- 1 STL6cO° TL tcl6co 1. “9c/6c0 2b Gyre ott 


9TLG6EO°T GUL6EO°T TeL6e0°T 9696E0°T ULLeco Tt 


LTLGEO*T OOL6GEO°T 9ICLEEO°T Enbecoo T 


OSS6e0°T veleco°T “tvloroe T 


O7PLEEO'T 


L80@Z2°0 = A = A Qe pEezeNnTeas ofp L(A) X] (A= A) 


07 squewtxorzdde sped wor pouTeqqo soTyIeT TeTubuetszzzedhy 


oe 
aya Fo ‘(~a) oO AQtTtqtadessns eyz Fo oepnattTdue sy fo sezyeUTIsY 


G°M ATIVE 


L 


oe ny 
is 


= fe 


’ oEre 


0.1 
Oeeeeo.£ SETCEC.L FBTEDOLS 


@iteco.L TitecoA ooresc.f sssee0.1 #2veeeis 
f - i 


. 


ggvéegis 2ftecd.c rsseeoSt sestedd 
DLTREQ.L BLTelOd  sSTREO 

eiveeo.. severg.£ 2eveeg.t 

pireco.t erteco.s Eveeeo.£ 


BETOES.L 


249 


6S8L°0 6 
z98L°0 «= 6S8L°O0-—sC«éTBL°O COB 

z98L°0 =6§- #98L°0 = 8S BLO l 

€98L°0 pl8L°0 T98L°0 9 

Z98L°0 «= «6 98L° 0 29840 S 
6S8L°0 ? 

6 8 L 9 S v 7\e 


£80Z22°0 = °A = A 3e pozenTerAe (LY/F-) OF SzueUTxordde 


aped woxrs peuTeqgo 8°T33eT zetnbuetazazzedhy ouq azoxy (On4/9F-) FO soqeutqasg” 


9° ATAVL 


cles aoe 


x 7 


--saer.t 0 epee. Oo  Spet.” 
€380.0..  -STST.0 dante 


ssat.o sasv.0  8z8T.0 
Saer.0.  . 86t.6 


250 


19 Cee 0 


Bcc -0 
87° 0 Lecce SoG 
9¢2°0 nec ECC 
LCG 0 0€c°0 EOC7O 
Gc2c-0 Svc 70 Sac O 
qec 0 VTc°O 


L80Z22°0 = Fam. A Ze pezentTease C/N O32 sjUueWTxordde sepeq wozy peuTeaAgqo 


SOTA3eT TeTNHuetzaqzzedAy suyA so MP ay/?n- ‘KBrSeusa TeOTATAS ey Fo soeqeutasyg 


£°M FTaVL 


ee ee Te 


eee : e - < - : .. ~ ) 
= le Rheba anigetagecne tart ber rep: or re 


Tac .0 
[e¢.0 


6g¢.0 


252 


TABLE K.8 


Estimates of the amplitude of the specific heat, Ai (Ko), 


of the hypertriangular lattice obtained from Pade approxi- 
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